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Abstract 

We consider an illiquid financial market where a risk averse investor has to liquidate a portfolio within a 
finite time horizon [0, T] and can trade continuously at a traditional exchange (the "primary venue") and 
in a dark pool. At the primary venue, trading yields a linear price impact. In the dark pool, no price 
impact costs arise but order execution is uncertain, modeled by a multi-dimensional Poisson process. We 
characterize the costs of trading by a linear-quadratic functional which incorporates both the price impact 
costs of trading at the primary exchange and the market risk of the position. The liquidation constraint 
implies a singularity of the value function of the resulting minimization problem at the terminal time T. Via 
the HJB equation and a quadratic ansatz, we obtain a candidate for the value function which is the limit 
of a sequence of solutions of initial value problems for a matrix differential equation. We show that this 
limit exists by using an appropriate matrix inequality and a comparison result for Riccati matrix equations. 
Additionally, we obtain upper and lower bounds of the solutions of the initial value problems, which allow 
us to prove a verification theorem. If a single asset position is to be liquidated, the investor slowly trades 
out of her position at the primary venue, with the remainder being placed in the dark pool at any point in 
time. For multi-asset liquidations this is generally not the case; it can, e.g., be optimal to oversize orders in 
the dark pool in order to turn a poorly balanced portfolio into a portfolio bearing less risk. 



1 Introduction 



In the last years, equity trading has been transformed by the advent of so called dark pools. These alterna- 
tive trading venues differ significantly from classical exchanges and have gained a considerable market share, 
especially in the US. Dark pools vary in a number of properties such as crossing procedure, ownership and 



accessibility (see Mittal 2008 and Degryse et al. 2009b for further details and a typology of dark pools). 
However, they generally share the following two stylized facts. First, the liquidity available in dark pools is 
not quoted, hence making trade execution uncertain and unpredictable. Second, dark pools do not determine 
prices. Instead, they monitor the prices determined by the classical exchanges and settle trades in the dark pool 
only if possible at these prices. Thus, trades in the dark pool have no or less price impact^] 

In this paper, we are concerned with the solution of a stochastic optimal control problem in continuous time 
arising in the context of optimal portfolio liquidation if an investor has access both to a classical exchange (also 
called "primary venue" or "primary exchange") and to a dark poolj^] We study a continuous-time model for 
optimal liquidation of a portfolio consisting of n assets within a finite time horizon [0, T] reflecting the trade-off 
between execution uncertainty of dark pool orders and price impact costs of trading at the primary venue. It 



complements the model by Kratz and Schoneborn |2012 , where (in particular) multi asset liquidation in dark 
pools is studied in discrete time; for single asset liquidation, the two models are connected via a convergence 
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[2003] and |Fong etaL] [2004] . 

he overall liquidity traded in dark pools in the US is strongly fragmented among approximately 40 different venues, see 



e.g., |Carrie| [2008] , Therefore, liquidity aggregation is a major issue. |Ganchev et al.| [2010] and |Laruelle et al.| |20lT] establish 
learning algorithms to achieve optimal order split between dark pools. Instead of analyzing the simultaneous use of several dark 
pools, we consider such an "aggregated" dark pool in our model which we call "the dark pool". 



result. The mathematical analysis of the continuous-time model is substantially different from the discrete- 
time model. In discrete time, the optimal liquidation strategy and its costs are given by a backward recursion 
obtained by standard dynamic programming; this recursion is rather unhandy and it is difficult to deduce 
properties of the optimal strategy from it. In continuous time, the liquidation constraint imposes a singular 
boundary condition for the value function and therefore the problem cannot be solved by standard methods 
of stochastic control. Instead, these methods had to be modified and extended. The optimal strategy and the 
value function are characterized via the "principal solution" of a matrix differential equation. We believe that 
the analysis of the differential equation and the resulting verification theorem for the optimization problem is 
mathematically interesting in its own right. Additionally, the differential equation is easier to analyze than 
the (discrete-time) backward recursion and we can deduce additional properties of the value function and the 
optimal strategy from it; in particular, we are able to show that for a two asset portfolio, the optimal strategy 
is monotone in the correlation of the assets. To our best knowledge, the mathematical framework established 
here and in Kratz and Schoneborn |2012 is the first within which optimal trade execution is analyzed for a 
multi asset portfolio using a classical exchange and a dark pool simultaneously. 

Our model for trading and price formation at the classical exchange is a linear price impact model. Trade 
execution can be enforced by selling aggressively, which however results in quadratic execution costs due to a 
stronger price impact. We model order execution in the dark pool by an n-dimensional Poisson process. Orders 
submitted to the dark pool are executed at the jump times of the respective components of the process. The 
split of orders between dark pool and exchange is thus driven by the trade-off between execution uncertainty 
and price impact costs. 

The investor aims to maximize her expected proceeds minus a quadratic risk term. For appropriate assump- 
tions on the price dynamics, this yields a linear-quadratic cost functional. However, the liquidation constraint 
implies a singularity of the value function at the terminal time T, which renders the solution of the optimiza- 
tion problem difficult. We approximate the constraint by a sequence of modified unconstrained optimization 
problems with increasing finite end-costs. Via a quadratic ansatz, the corresponding Hamilton-Jacobi-Bellman 
(HJB) equation suggests that the value functions of the modified optimization problems are quadratic forms 
for matrix-valued functions which are the solutions of initial value problems for a specific matrix differential 
equation. Explicit solutions for these initial value problems are not known. We establish a matrix inequality 
which allows us to apply a known comparison result for Riccati matrix differential equations in order to obtain 
closed form upper and lower bounds for the solutions of the initial value problems. The bounds are important 
for several reasons. Firstly, they enable us to prove a verification theorem for the modified unconstrained opti- 
mization problems with finite end-costs. Secondly, they imply the existence of the limit of the solutions of these 
optimization problems and thus yield a well-defined candidate for the value function of the original optimization 
problem with liquidation constraint. Thirdly, the limit of the bounds transfers to bounds for this candidate 
value function and we deduce a verification theorem for the solution of the optimal liquidation problem. 

Hence, the value function of the problem is a quadratic form for a matrix-valued function which is the 
"principal solution" of the matrix differential equation. The optimal strategy is linear in the portfolio position; 
it is determined by the value function and can be computed easily from a numerical solution of the "principal 
solution" of the matrix differential equation. This makes the model tractable for practicable applications and 
allows us to investigate the properties of the solution. When a position in a single asset is to be liquidated, 
the current asset position is at all times being offered in the dark pool, while it is liquidated in parallel at 
the primary exchange. The opportunity to trade in the dark pool leads to a slower liquidation at the primary 
exchange compared to a market without a dark pool. Traders hence need to fundamentally adjust their trade 
execution algorithms when a dark pool is introduced. It is not sufficient to use an algorithm that was optimal for 
execution at the exchange and to add a component that also places trades in the dark pool; instead, trading at 
the exchange needs to be adjusted during the entire trade execution time interval. While in a single asset setting 
the entire asset position is placed in the dark pool, this is not true if a multi asset portfolio is to be liquidated. 
If, e.g., the portfolio is balanced and thus only exposed to little market risk, then a complete liquidation of the 
position in one of the assets is unfavorable and thus only a fraction of the entire portfolio should be placed in 
the dark pool. This highlights again that overly simple adjustments to existing trade execution algorithms are 
exposed to potential pitfalls. For dark pools, the reluctance of traders to place balanced portfolios in a dark 
pool is an incentive to offer balanced executions in order to attract more liquidity. 

Our paper is connected to several lines of research. First, it builds on research on optimal trade execution 
strategies for a single trader in models with exogenously given liquidity effects. Several such models have been 



proposed for classical trading venues, e.g., Bertsimas and Lo 11998 , Almgren and Chriss 12001 , Almgren 12003 



Obizhaeva and Wang| |2006| , |Schied and Schoneborn| |2009| , |Alfonsi et al.| [2010] and |Schied et al.||2010| . We 



follow Almgren and Chriss|(200T] and assume a linear temporary price impact model for the primary venue. This 



choice yields a tractable model which nevertheless captures price impact effects. Linear price impact models have 
become the basis of several theoretical studies, e.g., 



Almgren and Lorenz 


|2007 , 


Carlin et al. 


|2007 , 


Schoneborn 



and Schied |2009 and Rogers and Singh 201?]]. Furthermore, they demonstrated reasonable properties in real 



world applications and serve as the basis of many optimal execution algorithms run by practitioners (see e.g., 
Kissell and Glantz 2003 ; Schack 2004 ; Abramowitz 2006 ; Leinweber 2007 ). The models above are limited 



to trading on one venue; they are concerned with the trade-off of execution costs versus market risk. By using 
a dark pool, the investor additionally faces order execution risk in our model. The trade-off of execution costs, 
market risk and execution risk is also apparent in the models o fjBayraktar and Ludkovski 2012 and Gueant 



et al. 2012 who consider optimal liquidation with limit orders^ In these models, the investor choses a limit 



price for her orders; the "probability" of execution is then dependent on this price and the left-over position can 
be liquidated at fixed costs at the exchange. In contrast to this, the dark pool orders in our model do not involve 
a limit price; they are executed at the price of the primary venue. The continuous-time model studied in our 
paper complements the discrete-time analysis of Kratz and Schoneborn 2012 for the specific case of linear price 



impact and no adverse selection. Klock et al. 2012 also work in a market model with exogenously specified 
liquidity characteristics for both a public exchange as well as a dark pool. Their focus is on the circumstances 
that can lead to price manipulation, while we focus on the quantitative and qualitative features of optimal trade 
execution strategies in markets without price manipulation. Therefore they investigate a more general class of 
price impact relationships, while our framework incorporates several aspects of market and investor behavior 
that shape optimal trade execution, such as risk aversion, multiple correlated assets and a dynamic usage of the 
dark pool. 

A second line of research connected to our paper focuses on the underlying mechanisms for illiquidity 
effects. Early examples include Kyle |1985 , Glosten and Milgrom [1985 and Grossman and Miller |1988| . In 
these models, price impact arises endogenously through the interplay of market participants. More recently, 
such models have been proposed t o analyze the competition between classical trading venues and dark pools. 
Hendershott and Mendelson 2000 analyze the interaction of dealer markets and a crossing networlj^jin a static 
one period framework where each investor buys or sells a single share. Their findings include that trading in a 
crossing network is cheaper than in dealer markets since the trader saves half the spread, but trade execution is 
uncertain in the crossing network. In a similar setting, Donges and Heinemann 2006 focus on game-theoretic 



refinements in order to remove the multiplicity of equilibria in Hendershott and Mendelson 2000 . Degryse et al. 



|2009a introduce a dynamic multi-period framework and compare the effect of different levels of transparency 
of the dark pool. In order to focus on the optimal order execution of an individual trader, we exogenously 
specify the liquidity properties of the market. Hence, all the models cited above, can shed more light on how 
the liquidity properties of the dark pool come about than our model can. On the other hand, by defining the 
model parameters exogenously, we can allow the investor to split her orders over time as well as between the 
two trading venues. While executing her order over time, she can dynamically react to the existence or absence 
of liquidity in the dark pool and adjust her trading strategy accordingly. Furthermore, we can take a multi 
asset perspective and investigate how the composition of basket portfolios influences optimal trade execution 
strategies. 

Methods of stochastic control is a third line of research that this paper is connected to. A standard reference 
for stochastic control with jumps is the book by 0ksendal and Sulem 2007 . The liquidation constraint yields 



a singularity of the value function at the terminal time T, and thus, the resulting stochastic control problems 
require non-standard considerations. For single asset liquidation and finite end-costs (hence no singularity 
of the value function at time T), the solution of the linear-quadratic control problem is well-known, see, e.g., 



0ksendal and Hu 2008 . Multi-dimensional linear-quadratic control with jumps is treated in the book by Hanson 
1 2007. However, our setting (even without the singularity of the value function) is not covered therein. The 
difficulty in our setting stems from the combination of the singularity of the value function and the fact that 
we consider multi-dimensional portfolios; thus the solution of the optimization problem involves the detailed 
analysis of a specific non-Riccati-type matrix differential equation, for which we establish existence results and 
upper and lower bounds of the solution by means of a novel matrix inequality. We are aware of only few other 

study optimal 



2010 



texts dealing with singular boundary constraints in liquidation problems. Schied et al 
liquidation without dark pools for CARA investors. In this case the optimal control problem does not include 
jumps. They carry out a verification argument for a candidate value function given in closed form. Naujokatj 



^In our model all orders in the primary venue are market orders whose execution is guaranteed; hence, the execution risk only 



applies to the dark pool orders. 

4 (Jrossmg networks are specific types ot dark pools that otter periodic matching ot buy and sell orders. 
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and Westray [2011 and Hoschler [2011 treat similar control problems with jumps. The focus of both texts is on 
trading with limit orders rather than with dark pools; they only treat single asset trading and obtain the single 
asset case of this paper as special cases of their respective settings. Kratz 2012 generalizes the single asset 



case of this paper by allowing for adverse selection; this renders the liquidation problem non-linear quadratic; 
he shows that in this case the value function is a "quasi-polynomial" of degree two instead of being quadratic; 
the solution for the multi-dimensional case it not known yet. 

The remainder of this paper is structured as follows. We specify the model setup and the optimization 
problem in Section [2] and obtain a candidate for the value function of the unconstrained optimization problem 
with finite end-costs via an initial value problem for a matrix differential equation. In Section [3] we state the 
main theoretical results of the article: the solution of the initial value problem (Section |3.1|, the solution of 



the unconstrained optimization problem (Section 3.2) and finally the solution of the constrained optimization 
problem (Section |3.3[ ) . We discuss the properties of the value function of the optimal strategy extensively for 
a single asset position and a portfolio consisting of two assets in Section |4j The main results of Section [3] are 
proven in Section [5] and the results of Section [4] are proven in Section [6] 



2 Model description 

For a fixed time interval [0, T], we consider the stochastic basis P, F = (^ 7 t)te[o,T]F| We investigate a 

market model where a risk averse trader with a personal risk aversion parameter a > has to liquidate a 
portfolio x G R n of n assets within a finite trading horizon [0, T]. The investor has the possibility to trade 
simultaneously at a traditional exchange and in a dark pool, facing the trade-off of paying market impact costs 
in the traditional exchange against uncertain execution in the dark pool. 

We specify the price dynamics at the primary exchange in Section |2.1| and the dynamics of order execution 
in the dark pool in Section |2~2| In Section [273] we define admissibility of trading strategies. In Section [2~4] we 
specify the trading objective and show that the resulting optimization problem is equivalent to a linear-quadratic 
optimization problem. Heuristic arguments suggest that the value function of the optimization problem is 
singular at terminal time T because of the liquidation constraint. Hence, we introduce a modified optimization 
problem where we drop the liquidation constraint and approximate it by finite end-costs for a portfolio not 
liquidated by time T as an intermediate step. In Section [2. 5 1 we derive a candidate for the value function of the 
modified optimization problem via a quadratic ansatz and the corresponding HJB Equation. 

2.1 The primary exchange 

In absence of transactions of the investor, the fundamental asset price at the primary exchange is given by an 
n-dimensional stochastic process P. 

Assumption 2.1. We assume that P satisfies the following properties, 
(i) P is a square-integrable cadlag martingale. 

(ii) The covariance matrix of P is constant in time, i.e., for all t G [0, T], i,j = l,...,n, 

Cov(Pi(t),Pj(t)) =iEij 

and E = (Sij)ij=l,...,n € R nxn - 

As the time-horizon for portfolio liquidation is usually short, i.e., several hours or a few days, the martingale 
property in Assumption 2.1 (i) does not constitute a major divergence from reality^] Once the trader becomes 
active on the primary exchange, she influences the market price P. We assume that the trader at the primary 
exchange can only execute trades continuously, i.e. that her trading activity on the primary exchange is abso- 
lutely continuous and can hence be described by her trading intensity £(t) with t G [0,T). This trading in the 
traditional exchange generates price impact, which we assume to be temporary and linear in the trading rate 
£(t). Given a strategy (£(£))te[o,T]i the transaction price at time t G [0, T] is given by 

P(i)=P(*)-A£(f), 

5 The filtration is generated by the involved random processes and is specified at the end of Section 2.2 

6 We do not model the period in which the investor has held the assets before time zero. For this period, the martingale 
assumption is naturally a worse approximation of reality. 
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where A £ R" x ™ is a positive definite matrix constant in time. Execution of the trades at the primary exchange 
is certain; we hence consider only market orders and no limit orders. 

By assuming linear price impact for the primary venue, we follow Almgren and Chriss |2001 . This choice 
yields a tractable model which nevertheless captures price impact effects. Linear price impact models have 
become the basis of several theoretical studies, e.g., 



Almgren and Lorenz 


|2007 , 


Carlin et al. 


|2007 , 


Schoneborn 



and Schied 2009 and Rogers and Singh 2010 . Furthermore, they demonstrated reasonable properties in real 
world applications and serve as the basis of many optimal execution algorithms run by practitioners (see e.g., 



Kissell and Glantz 2003 ; Schack 2004 ; Abramowitz 2006 ; Leinweber 2007|). In contrast to Almgren and 



Chriss |2001 , the price impact in our model is purely temporary; |Klock et al.| |2012 analyze the influence of 



permanent price impact on the existence of market manipulation strategies. Such an analysis is not in the scope 
of our article; hence we allow only for temporary price impact Q 

2.2 Order transaction in the dark pool 

In addition to the primary exchange, the trader can also use a dark pool. Dark pools often have rather complex 
order allocation mechanisms; most of them use some sort of a pro-rata or time-priority rule for matching orders 
from opposite sides of the market. Here, we consider a dark pool with a time-priority matching rule: the 
investor's order rj enters a queue and is matched with liquidity from the opposite side of the market (if there is 
any) once it has reached the front of the queue. 

We allow for continuous updating of the orders rj(t) in the dark pool at any time t £ [0, T]. Orders for the 
z th asset in the dark pool are executed fully at the jump times of the i th component of an n-dimensional Poisson 
process 

7r = (7Ti, . . . , 7r„) with intensities 6i, . . . , 9 n > 0, respectively. 

Else, the orders are not executed at all. This mechanism implies two main simplifications of reality which 
allow a thorough mathematical analysis of the model. First, we exclude partial execution; the probability of 
execution does not depend on the size of the order. Second, we assume independence of the increments of the 
dark pool liquidity. On the other hand, the resulting model captures the stylized facts outlined above and in 
the introduction; we believe that it constitutes a sufficiently well approximation of reality for our purposes. 

Assumption 2.2. We assume that tt satisfies the following conditions. 

(i) 7Ti, . . . , 7r„ are independent. 

(ii) 7r and P are independent. 



By Assumption 2.2 (ii) we exclude correlations between dark pool liquidity and the fundamental asset price 
in the market. In particular, this rules out adverse selection^] For single asset liquidation, adverse selection 
was incorporated by Kratz |2012 . This renders the liquidation problem non-linear quadratic; he shows that in 
this case the value function is a "quasi-polynomial" of degree two instead of being quadratic. A solution for the 
multi-dimensional case it not known yet. 

While the dark pool has no impact on prices at the primary venue, it is less clear to which extent the price 
impact of the primary venue A£(t) is reflected in the trade price of the dark pool. If for example the price 
impact is realized predominantly in the form of a widening spread, then the impact on dark pools that monitor 
the mid quote can be much smaller than A£(i). We will make the simplifying assumption that trades in the 
dark pool are not influenced by the price impact at all, i.e., that they are executed at the fundamental price 
P. If alternatively the transaction price in the dark pool is the price P at the primary exchange including the 
trader's price impact, market manipulation strategies can become profitable unless the parameters are chosen 
with great care, as has been shown by Kratz and Schoneborn 2012 for the discrete-time case. For a detailed 



discussion see also |Klock et al. 
dark pools. 



2012 



who analyze the circumstances that can lead to price manipulation in 
We are now ready to specify the filtration (J-t)t as the completion of ( cr(P(s), 7r(r)|0 < s <t, < r < t)) t . 



7 As the price is not influenced permanently by the investor, the term price impact might be misleading. Alternatively, we could 
rename the quadratic costs caused by trading in the primary venue as the execution costs of the investor; indeed, the costs of 
trading can be interpreted to include such different effects as price impact and transaction costs. However, the term execution 
costs does not capture the possible impact of trading in one asset on the price of another asset which we include by allowing the 
off-diagonal elements of A to be non-zero (cf. the discussion of cross price impact in Section |4.2.2| . 

8 Here, adverse selection refers to the phenomenon that liquidity seeking traders find that their trades in the dark pool are usually 
executed just before a favorable price move, i.e., exactly when they do not want them to be executed since they miss out on the 
price improvement. 
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2.3 Admissible trading strategies 

Let t £ [0, T) be a given point in time and x £ R" be the portfolio position of the trader at time t. The trader 
has the possibility to trade asset k in the traditional exchange with trading intensity £fc(s) at time s £ [t,T) 
and to place orders r)k{s) in the dark pool at time s. 
We call a 2n - dimensional stochastic process 

(u(s)) se[ uT) = (f(s),»7(s)) se [ t) T) 

a trading strategy if £ is progressively measurable and r) is predictable. Given a trading strategy u, the portfolio 
position at time s £ [t, T) is given by the following controlled stochastic differential equation: 

dX u ( S ) = -t(s)ds~r]( S )d7r(s) 
X u (t) = x U 

such that the left hand side in is well-defined. 

For technical reasons we require all trading strategies to fulfill the following conditions. 

Definition 2.3. Let t £ [0, T) and x £ R™ be fixed. Let u = (v,(s)) a £\t,T) = ((£( s )i v( s )))sg[o,t) be a trading 

strategy, i.e., £ is progressively measurable and rj is predictable. 

(a) We call u an admissible trading strategy if it fulfills the following conditions. 

(i) The Stochastic Differential Equation possesses a unique solution on [t,T). 

Y llfOOIftfe] <oo, E[J \\v(s)\\ S 2 ds 



E 



< 00. 



(iii) If Oi = 0, then ^(s) = for all s £ [t,T). 

We denote the set of admissible trading strategies by A(t). 

(b) We call u £ A(t) an admissible liquidation strategy or just liquidation strategy if additionally 

(iv) lim s _j.T- X u (s) = a.s. 

and denote the set of admissible liquidation strategies by A(t,x). 



Let us shortly comment on Definition 2.3 Condition (ii) is required for the moment bounds in Lemma 5.5 
which is in turn needed for the verification later. Condition (iii) is needed in order to ensure uniqueness of 
optimal trading strategies: if 9i — 0, no additional gain can be achieved by non-zero dark pool orders. If the 
portfolio is liquidated at constant speed at the primary exchange only, i.e., £(s) = jttj, tj(s) — for s £ [t,T), 
Definition |2.3| (in particular the liquidation constraint (iv)) is satisfied and hence A(t,x) ^ 0. 



Remark 2.4. We expect that the stochastic control problems we solve are such that the optimal control is of 
Markovian form (see, e.g., the book by \0ksenda\ \2007l, Theorem 11.2.3): 

u{s) = (£(*),»?(*)) = (S(8,X(s))M*,X(*-)) 

for deterministic functions £,r) : [t,T) x R n — > R™. The deterministic initial value problem 

X' = -a-,X), X(t)=x (2) 

possesses a unique solution on [t, T) if ||£(s, y)\\ < /(s)||y|| + g(s) on [t, T) x R n for f,g£ C([t, T)) and £(s, •) 
is locally Lipschitz (e.g., C 1 )f\ Let £ : [t, T) x R" -> R" fulfill these conditions and let -q : [t,T) x R n -> R". 
We can pathwise construct the solution of the Stochastic Differential Equation ([TJ inductively by interlacing the 
jumps (see, e.g., App lebaum\ Example 1.3.13): as the n Poisson processes are independent, they jump at 



distinct times almost surely. Let (Tj)j>i be the jump times of it such that t =: r < T\ < . . . almost surely, and 
let X be the solution of the Initial Value Problem Q on [T,*,Ti+i A T) with initial value x — X{t,i) for i g I 
such that Ti < T. For rj+i < T and A7T/ c (ri + i) > 0, we set 

x {n+i) ■= X(r i+ i-) - r)k(T i+ i,X(Ti+i-))e k , 

where is the k th unit vector. 



3 This follows, e.g., by Peano's existence theorem and Gronwall's inequality. 
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2.4 Cost functional 



The proceeds of selling the portfolio x e R™ during [t,T] according to the strategy (u(s)) s = (£,(s),rj(s)) s E 
A(t, x) are given by 

<Kt,x,u):= j £( S ) T (P( S )-A£( S ))d S + j v( S ) T P(s)dn(s). 

The first term in the above equation represents the proceeds of selling at the primary exchange at a price of 
P(t) = P(s) — A£(s), while the second term accounts for the proceeds of selling in the dark pool at the un affec ted 
price P(s). Applying integration by parts and using X u (s) = x — £,(r)dr — r](r)dTr(r), Assumption 
the fact that P and it are independent (Assumption 2.2 (ii)) and the liquidation constraint (Definition 
we obtain 



2.1 



2.3 



(i), 



4>(t, x, u) 



£(s) T A£(s)ds + x T P{t) + [ X u (s~)dP{s). 
t Jt 



This yields (cf. Assumption 2.1) 



E [(f)(t, x 7 uj] = x T P(t) — E j £(s) T A£_(s)ds 

Instead of maximizing expected proceeds, we can thus equivalently minimize expected price impact costs. We 
assume that the trader is not only interested in expected liquidation proceeds, but in addition also wants to 
minimize risk during liquidation. We incorporate both aspects in the following cost functional p^] 



J(t,x,u) :=x T P(t) -E[(j)(t,x,u)] +E a / X u (s) T Y,X u (s)ds 



= E 









f(as),x u ( s ))ds 



where / : R™ x R™ — > R is given by /(£, x) :— £ T A£ + ax T Hx . The first two terms in the cost functional capture 
the expected liquidation shortfall, while the last term is an additive penalty function a X u (s) T Y,X u (s)ds 
which reflects the market risk of the portfolio (recall that a > is the risk aversion parameter of the investor) ; 
it penalizes slow liquidation and poorly balanced portfolios. It does not incorporate execution risk; hence the 
investor is only risk averse wit h respect to market risk but not with respect to execution risk. In discrete time, 
Kratz and Schoneborn 2012 argue that for realistic parameters, market risk outweighs execution riskj^J we 
hence expect that the inclusion of risk aversion with respect to execution risk would not change the optimal 
strategy in essence. For deterministic liquidation strategies without dark pools, the risk term reflects the variance 



of the liquidation costs (see Almgren and Chriss 2001 ). In this case, minimizing a mean-variance functional 
of the liquidation costs over all deterministic strategies is equivalent to maximizing the expected utility of the 
proceeds of an investor with CARA preferences over all strategies (see Schied et al. |2010| ). 

We assume that the trader aims to minimize the cost functional and hence considers the following optimiza- 
tion problem PI 

v(t,x):= inf J(t,x,u). (OPT) 

«eA(t,x) 

Note that the optimization problem is well-defined and the value function satisfies v(t,x) < oo for t < T 
(consider, e.g., constant liquidation exclusively in the primary exchange). Because of the liquidation constraint 
(cf. Definition 2.3 (iv)), we expect the value function to fulfill 



lim v(s.x) 







if x = 
else, 



10 Both mean and variance of execution costs are often used as measures of execution quality. The cost functional J is inspired 
by such mean variance measures. An alternative approach is the maximization of expected utility: J(t,x,u) = E[U(4>(t, x, u))] for 
some utility function U. In this alternative set-up analytical solutions are unfortunately not directly available through the methods 
presented in this paper and are hence left for f uture research. 

11 The illustration and argument provided in |Kratz and Schoneborn] |2012] applies for the continuous-time case in the same way 
as it applies to the discrete-time setting in which it is presented. 

12 An alternative interesting set-up is to consider trade execution under minimum proceeds constraints. Such problems have 
recently been addressed using theory about stochastic target problems (see, e.g., |Bouchard et al. 2009 and Bouchard and Dang 
|2012| ). 
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i.e., v has a singularity at the terminal time T. Because of this singularity, non-standard considerations are 
necessary for solving the Optimization Problem (OPT I via a verification argument using the HJB equation. 



As an intermediate step, we hence weaken the liquidation constraint by allowing for all strategies u G 
and by penalizing non-liquidation by finite end-costs. More precisely, for / > and u — (£,rj) € A(f), we define 
the following cost functional 



J(l, t, x, u) := E 
The resulting optimization problem is 



/(£ (s), X u (s))ds + I ■ X U {T) T IX U {T) 



v(l,t,x) := inf J(l,t,x,u). 
ueA(t) 



(OPT) 



The Optimization Problem (OPT I mainly serves as an approximation of the optimization Problem (OPT I. 
However, it is also of interest itself: the penalization term I ■ X M (T) T I X u (T) can be considered as the liquidation 
cost of the left-over position. Note that the identity matrix in the term can be replaced by any positive definite 
matrix reflecting this interpretation (e.g., A) without changing any of the proof s significantly. 

In the following, we solve the unconstrained O ptimiza tion Problem (OPT) first (Section 3.2 1. Then, we 
show that the solution of the Optimization Problem ( OPT I converges to the solution of the original constrained 



Optimization Problem ( OPT ) as I — > oo (Section |3.3 1 . 
2.5 Hamilton-Jacobi-Bellman equation 



In this section we derive a candidate for the value function of the Optimization Problem ( OPT I . Heuristic 
considerations suggest that it should satisfy the following HJB equation (see, e.g., Oksendal and Sulem |2007| ): 



dw 



Of 

w(T, x) — Ix 1 x. 



(t,x) = sup di(w(t, x) - w(t, x - T] T ei))+\7 x w(t, x)£-f(£,x) 

u=K,i,)£R"xE" L i=1 
T 



(HJB) 



The linear-quadratic form of the cost functional suggests that the value function is quadratic. Assuming that 
the above guesses are correct, the following proposition provides candidates both for the value function and for 
the optimal strategy. 

Proposition 2.5. Let I > and assume that the initial value problem for a matrix differential equation 

C = C T A- 1 C + C T CC-aE, C(T) = II, 

where C(l,t) := diag ( — s ^ possesses a positive definite solution C(l,t) — (cij(i, £))i,j=i,...,n on [0, T]. Then 

w(l, t, x) := x T C(l, t)x 



satisfies the HJB Equation (HJB I with maximizer u* = (£*,?7*) for 

C :=£*(l,t,x) := A^CXMK V* := »7*(l,t,ar) := C{l,t)C(l,t)x, 



where C(l, t) := diag ( ■ 



If there exist ix, . . . , if, G {1, . . . , n} such that 9i j =0 (j = 1, . . . , k), then Tjj. can be chosen arbitrarily. Up 
to arbitrary choices of r}i . , the maximizer is unique. 

Proof. Th e asse rtion follows directly from plugging the quadratic ansatz w(l,t,x) = x T C(l,t)x into the HJB 
Equation (HJB I: the resulting function can be maximized by standard calculus. □ 
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3 Main results 



Proposition 2.5 suggests that the solution of the Optimization Problem (OPT I solves the initial value problem 



for the matrix differential equation 



where 



C = C T A _1 C- 
C{T) = II, 

C := diai 



C 1 CC - aS 



(— )• 



(3) 



(4) 



In the remainder of the section, we state the main results of the article. In Section 3.1 we show that ([3| admits 
a unique solution C on [0,T] and establish approp riate upper and lower bounds for C. Subsequently, we deduce 
the solution of the Optimization Problem (OPT I in Section 3.2 and as a limit of this (as I — > oo) the solution 



of the Optimization Problem (OPT I in Section 3.3 Proofs of these results are presented in Section [5] 
Before we proceed, we introduce the following notations. 

Notation 3.1. (i) For symmetric matrices A and B we say A > B (A > B) if A~B is positive (nonnegative) 
definite. 



(ii) We denote the smallest and the largest eigenvalues of a real- symmetric matrix A by a n 
respectively. 



and 



(Hi) We define the following nonnegative definite matrix: D := \/A l£\/A . 

3.1 Solution of the Initial Value Problem ^ 

It is not immediately clear that the Initial Value Problem ([3| possesses a positive definite solution on the whole 
interval [0, T] for n > 2. For n — 1, it reduces to C = ^- + Q\C - a£, C(T) = I. This is an initial value 
problem for a scalar Riccati differential equation with constant coefficients, whose unique solution is explicitly 



known and exists on the whole interval [0, T] (cf. Section 4.1 1. For n > 2, the following theorem establishes the 
existence and uniqueness of the solution of (|3| . 



Theorem 3.2. Let 8i > for i = 1, . . . , n, 6 = @i an d I > where 

»=i 



; := max|A max ^y^ 



i(\/admax)| 



(5) 



(cf. Notation 3.1). Then the Initial Value Problem ^ possesses a unique solution C{1,-) on {— oo, T]. The 
solution is symmetric for all t G (— oo, T] and 



< P(l, t) < VA^C(l, t)VA^ < Q(l, t), 
where P and Q are the solutions of the initial value problems 



pi = p2 



6P 



P{T) 



/ and Q' = Q 2 



Q(T) 



respectively. 



(6) 



Remark 3.3. (i) The solutions of the initial value problems for Riccati matrix differential equations in ^ 
exist on the whole interval [0, T] and can be computed in closed form (cf. Equations (211 - (241/ For 
technical reasons, we prefer to establish bounds for \/A _1 C\/A _1 instead of bounds for C. P and Q are 
constructed in terms of multiples of the identity matrix and hence commute with all matrices. Therefore, 
they transfer to bounds of C directly by multiplying them with A. 



(ii) The bounds of C are an essential component for the proof of T heorem 3.2 Additio nally , they are req uired 

and Theorem 



for all key steps of the solution of the Optimization Problem (OPTl (Proposition 3.4 



and of the solution of the Optimization Problem (OPTl (Theorem \3.6\ Theorem 3.8 and Theorem 3.9) 



3.5) 
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3.2 Solution of the Optimization Problem (OPT) 



Combining Proposition |2.5| and Theorem |3.2| we obtain well-defined candidates b oth for the value function 
(x T C(l, t)x) and for the optimal strategy u* = (£*,7?*) of the Optimization Problem (OPTl. The latter is given 

by 



n*{l) := T)*(l,t,x) := IC(l,t)C(l,t)x, 
where / = (ei,j)i,j=i,,,., n is the diagonal matrix with 



(7) 
(8) 



ei-, 



if Oi > 
else 



and 



C(l,t) := dia; 



iag (c, i(Z,i)) 



The following Proposition confirms that u* is admissible. 

Proposition 3.4. Let I > l forl as in Equation (J5) and (t, x) G [0,T)xR". Thenu*(l) = (£*(0>»7*(0) € A(t), 
where £*(Z) and T)*(l) are as in Equations Q and ([8j ; respectively. 

Finally, we obtain the solution of the Optimization Problem (OPT). 

Theorem 3.5. Let I > Iq for Iq as in Equation (|5| and let C(l,t) be the unique solution of the Initial Value 
Problem pj. Then the value function of the Optimization Problem (OPTl is given by 

v(l, t, x) = x T C(l, t)x 

and the P (§5 A - almost surely unique optimal strategy is given by u*(l) as in Equations ^ and (pj. 



3.3 Solution of the Optimization Problem (OPT) 

Intuitively infinite end-costs should force the controlled process X u (s) to approach zero as s — > T-; furthermore 
the soluti on of the Optimization Problem ( OPT I should be the limit of the solution of the Optimization 
Problem (OPTl as / — > oo. The following theorem confirms (in particular) that this limit is well-defined. 

Theorem 3.6. Let t e [0,T). 

(i) The element-wise limit of the value function matrix 

C(t) := lim C(l,t) 

l— >oo 

exists on [0,T), and C(/, •) converges compactly to C on [0,T). Furthermore, lim/_ i . 00 c m i n (Z, T) = oo. 

(ii) C solves the matrix differential equation 

C = C T A- 1 C + C T CC-aY, (9) 

on [0, T) with boundary condition lim^j^ c m in(s) — oo. Moreover, the following inequalities hold. 

< P(t) < ^fk^Cit)*/! 1 ^ < Q(t), (10) 

where P(t) := lim^oo P(l, t) and Q(t) := lim/^oo Q{1, t). 

Remark 3.7. For Riccati matrix differential equations, there exists a unique solution F with lim s _j.y_ /min( s ) 
X This solution is called the principal solution (see, e.g., Coppel \1911\ ). In this spirit, C is the principal 
solution of the Matrix Differential Equation (J9|) . Note however that it is not entirely clear that C is the only 
solution of (J9j) satisfying lim s ^T- c m in(s) = 00 since Q is not a Riccati matrix differential equation. 

By Theorem |3.6[ we also obtain the existence of the limits of the optimal strategy: 

:=£*(*, a) := lim C(l,*,x) = A^Ctyx, rf := rf(t,x) := lim rf{l,t,x) = IC(t)C(t)x. (11) 



It turns out that u* := (£*, rf) is an admissible liquidation strategy, in particular that it satisfies the liquidation 
constraint 

lim r(s) = for X*(s) := X u * (s). 
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Theorem 3.8. Let t £ [0,T), x £ R n and u* = (£*,»?*) /or £* and rf as m (TTTJ . Then u* £ A(t,x) 



We are now ready to present the main result of this article: the solution of the Optimization Problem (OPT I 



Theorem 3.9. The value function of the Optimization Problem (OPT I is given by 

v(t,x) = x T C(t)x 

for all t £ [0,T), x £ R™ and 

JO ifx = i) 
km v{s,x) = < 
s-yr- oo else. 



The 



X - almost surely unique optimal strategy is given by u* — (£*,r)*) as in (111 



4 Properties of the value function and the optimal strategy 



The characterization of the solution of the Optimization Problem ( |OPT[ ) enables us to analyze the properties 
of the optimal strategy and the value function in detail^] For single asset liquidation (n = 1) the Initial Value 
Problem (J3j can be solved in closed form as the differential equation is a scalar Riccati equation with constant 
coefficients. This allows us to prove monotonicity properties of the value function and the optimal strategy in 
Section [4~T] In Section [472] we discuss the multi asset case by analyzing a portfolio of two assets. Although a 
closed form solution of the Initial Value Problem ^ is not known in general for n > 2, it is possible to derive 
analytical results about the dependence of the value function and the optimal strategy on the model parameters, 
in particular on the correlation of the assets. We illustrate that overly simple adjustments of existing trading 
algorithms for optimal liquidation without dark pools can have undesirable properties. The proofs of the results 
of this section are presented in Section [6] 



4.1 Single asset liquidation 

We let n — 1 and set 9 = Oi . The solution of the Initial Value Problem (|3| is given by 



C(I,t) = T coth(-(T -*)+«(!)) 



A9 



where 



n(l) := arcoth 



lri.S 

A 



for 9 > or aY, > and C(l,t) — v, a for 9 = aE = 0. In order to highlight the dependence of the value 
function on the parameters 9 7 A and aS, we define, e.g., 



C{t;9) := C(t) 



lim C(l,t) 

I— too 




(T-t) 



AO 
2 



if 9 > or aS > 
if 9 = aT, = 0, 



in particular 



C(i;0) 



VaEAcoth (y/^(T - t)\ 



A 
T-t 



if aS > 
if a£ = 0. 



We will apply similar notations throughout Section [4] to make the dependence of other model components (such 
as the optimal strategy) on the respective parameters explicit whenever this clarifies the exposition. For 9 = 0, 
we obtain the special case of optimal liquidation without dark pool, see Almgren and Chriss| |2001| for the 



discrete-time case and Schied et al.||2010] for the continuous-time version. 



No transaction costs must be paid in the dark pool; intuitively, the investor should hence try to liquidate as 
much as possible in the dark pool. Indeed, we have 

r)*(t,x;9) = C(t;9)C(t;9)x = x 



13 We limit the analysis to the Optimization Problem ( |OPT| l. Most of the properties transfer directly to the Optimization 
Problem I OPT i with the same or similar proofs. 
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for 8 > by Theorem 3.9 i.e., it is optimal to always place the full remaining asset position in the dark poolp*] 



note that the execution of the dark pool order immediately stops the trading activity by linearity of the optimal 
strategy in the position. 



In Section 4.1.1 we discuss the dependence of the optimal strategy and the value function on 8. Subsequently 



we analyze the dependence on the price impact A and the risk parameter aS in Section |4.1.2| 
4.1.1 Dependence on 8 

We expect it to be optimal to slow down trading in the primary venue initially as the trader hopes to trade 



cheaper in the dark pool. This intuition is confirmed by Proposition 4.1 (ii) and (iii) below. In order to state the 
property rigorously we first denote the optimal trading trajectory until execution in the dark pool by X(-;8), 
i.e., X is the solution of the linear initial value problem X' — —£*(•, X; 9), X(0) = x, where x is the initial asset 
position at time zero. Then for t G [0,T), 

9) = xexp ( - /' gfeS*) = ^ - '») 7 «') x . (!2) 
v Jo A ' sinh(fT) 

We obtain the following monotonicity properties. For simplicity of exposition we assume a£ > 0. Similar 
results hold for the simpler case aS = (cf. also the right hand graph of Figure [T]) . 

Proposition 4.1. (i) Fort G [0,T) ; C(t;9) is strictly decreasing in 9. 

(ii) For x > and t £ (0, T], X]8) is strictly decreasing in 9. 

(iii) For x > and t £ (0,T), X(t;9) is strictly increasing in 9. 

(iv) For x > and t G (0,T), the expected asset position if the optimal strategy is applied, E[X*(t;8)], is 
strictly decreasing in 9. 

(v) For aY, > 0, the risk costs aE • E [ X* (t; 9) 2 dt\ are strictly decreasing in 9. 

(vi) The impact costs A • E[ f Q T £* (t. X*(t; 9); 9) 2 dtj are strictly decreasing in 8. 

Let us shortly comment on Proposition |4.1| The fact that the overall costs are decreasing in 9 (part (i)) is 



quite intuitive and can be deduced directly from the definition of the cost functional J. Proposition 4.1 (iv) 



states that the introduction of a dark pool decreases the expected asset position despite slower initial trading in 



the primary exchange (parts (ii) and (iii)). Parts (v) and (vi) of Proposition 4.1 confirm that the introduction of 
the dark pool decreases both the impact costs component and the risk costs component of the value function^ 
We illustrate these properties in Figure [l] In the left picture we consider a risk neutral investor and in the 
right picture a risk averse investor. The optimal trading trajectories for trading with dark pool are displayed 
by the thick solid lines. In the displayed scenario the dark pool order is executed at time r. After execution in 
the dark pool, the liquidation task is finished. The dotted lines denote the scenario where the dark pool order 
is not executed during the entire trading horizon and the thin solid lines refer to the expected asset positions. 
We contrast the optimal strategy with dark pool by the optimal strategy without dark pool (dashed lines). As 
shown in Proposition |4.1| the dark pool slows down trading in the primary venue initially. Nevertheless, the 
expected position is smaller than the trading trajectory without dark pool. 

4.1.2 Dependence on A and aE 

It follows directly from the definition of the cost functional J that the value function is strictly increasing both 
in the impact costs parameter A and in the risk costs parameter aE. For the optimal strategy, impact costs 
and risk costs have conflicting influences: while larger impact costs yield a reduction of the trading intensity, 
larger risk costs yield faster trading (cf. also the difference of the left and the right picture of Figure 0. We 
summarize these findings in the following proposition. 



14 We want to remark that this property is sensitive to the assumption that P is a martingale (Assumption 2.1 (i)). If, e.g., the 
investor holds a long position in the asset and P has a positive drift, she should be reluctant to sell her entire position too early. 
Therefore, we expect that a drift changes this property. Similarly, if adverse selection is included into the model, the property 
r)*(t,x) = x does not hold; this was shown bylKratzl [20121. 



It is an interesting effect that this is not necessarily the case in the discrete-time setting of Kratz and Schoneborn| |2012| where 



the risk costs can be increasing for small B. This effect is lost by letting the length of the trading periods tend to zero (see |Kratz| 



poTT ) 
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Size of asset position Size of asset position 




Figure 1: Optimal trading trajectories for a risk neutral (left picture) respectively a risk averse trader (right 
picture). The thick solid lines denote the optimal trading trajectory with dark pool in the scenario where the 
dark pool order is executed at time r. The dotted lines denote the scenario where the dark pool order is not 
executed during the entire trading horizon and the thin solid lines refer to the expected asset positions. The 
dashed lines denote optimal liquidation without dark pools. i = l,T = l,A = l, 9 = 4 and for the right 
picture a = 6, E = 1. 



Proposition 4.2. (i) Fort G [0,T), C(t; A, aE) is strictly increasing in A and in aE. 

(ii) Let t G [0,T) and x > be fixed. Then £*(t, x; A) is strictly decreasing in A. Consequently, X(t;A) is 
strictly increasing in A for t G (0,T). 

(Hi) Let t G [0,T) and x > be fixed. Then £* (t, x; aE) is strictly increasing in aE. Consequently, X{t;dS) 
is strictly decreasing in aE for t G (0, T). 



4.2 Portfolio liquidation 

If a risk averse investor has to liquidate a portfolio of multiple assets (n > 2) , then correlation between the assets 
comes into play. Depending on whether the portfolio is well (poorly) balanced, it is intuitively optimal to place 
orders in the dark pool which are smaller (larger) than the remaining portfolio position. In the first case, the 
risk costs of the portfolio are small and therefore, the trader will not risk to lose her balanced position by the full 
execution of her dark pool order for only one of the assets; hence her orders are smaller than the remainder of 
the position. In the second case, it might be optimal to place orders in the dark pool which are larger than the 
remainder of the position for risk mitigation reasons. This illustration suggests that overly simple adjustments 
of existing trading algorithms for optimal liquidation without dark pools can have undesirable results. 



In Section 4.2.1 we verify the above intuition by analytical results about the dependence of the value function 
and the optimal strategy on the correlation of a portfolio of two assets; we then deduce the general structure 
of the optimal strategy dependent on whether the portfolio is well or poorly balanced. We also discuss the 
dependence of the value function and the optimal strategy on the price impact parameter and the execution 
intensities in Section |"4.2.2| and |4.2.3| respectively. 

As a prerequisite, we introduce a characterization of the optimal dark pool order which exploits that the 
jump times of the Poisson processes are almost surely distinct. Intuitively, an execution of the optimal dark 
pool order for asset i should bring the position in asset i to its optimal value given unchanged positions in all 
other assets j 7^ i. The following Proposition confirms this intuition (see also Naujokat and Westray [2011 for 



a similar result). 

Proposition 4.3. Let t G [0, T], x G R n be the portfolio position at time t and i = 1, ...,n. Then, 

v(t, x — rj*(t, x)ei) = min v(i, x — rjei). 

?7SR 
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4.2.1 Dependence on correlation 



We will see that dark pool trading is sensitive to the correlation of price increments. In the following, we discuss 
the case n = 2, p*| In order to simplify the exposition, we assume that there is no cross asset price impact 

Ai \ _ / a\ po- x a 2 \ 



a= r: 1 , £= ^ *~ A a . (13) 

For the purposes of this section, we assume that the variances <j\ and a 2 of the two assets as well as the risk 
aversion parameter a are strictly positive. 

If the correlation of the two assets is positive (p > 0) , a portfolio consisting of a long position in one asset 
and a short position in the other asset is more desirable than long positions (or short positions) in both assets; 
in the former case, a part of the risk of each asset is hedged by the other asset. Conversely, if p < 0, it is more 
desirable to have long (or short) positions in both assets. 

Definition 4.4. A portfolio x = (xi, x 2 ) T (xi,x 2 ^ 0) is well diversified if either the signs of the positions are 
equal (sgn(xi) = sgn^)) and p < or if the signs of the positions are different and p > 0. Otherwise, the 
portfolio is poorly diversified. 

Proposition 4.5. Let t G [0, T) and x±,x 2 7^ 0. Then v(t,(xi,x 2 ) T ) < v(t, (x\, — x 2 ) T ) if and only if the 
portfolio x is well diversified; v(t, (xi, x 2 ) T ) > v(t, [x\, — x 2 ) T ) if and only if the portfolio x is poorly diversified. 

We can further specify the dependence of the value function on the correlation if the portfolio is well 
diversified. 

Proposition 4.6. Let t G [0, T) and x be well diversified. Then v(t,x) is strictly decreasing in \p\. 

The left picture of Figure [2] illustrates the dependence of the value function on the correlation p for a portfolio 
that is long in both assets. For p < 0, this portfolio is well diversified and the value function is increasing in 



p (i.e., decreasing in \p\) in line with Proposition 4.6 For p > 0, the portfolio is poorly diversified. This leads 
to elevated liquidation costs for small positive p. For large positive p, the increased risk costs of the current 
portfolio are outweighed by the (projected) smaller risk costs of a future well diversified portfolio (e.g., after the 
execution of an order in the dark pool; cf . also Proposition |4.8| (ii) below) . The opportunity of risk reduction 
results in a decrease of the value function for large values of p in the displayed example. 




Figure 2: Dependence of the value function (left picture), the optimal trading intensity for the second asset 
(middle picture) and the optimal dark pool order for the second asset (right picture) on the correlation p. 
x\ = x 2 — 1, T — 1, Ai = 3, \ 2 — 0.2, 6i = 0.5, 2 — 5, a — 4 and o\ = o 2 = 1. 

We have the following symmetry and monotonicity properties for the entries of the value function matrix C. 

Proposition 4.7. Let t G [0, T) and denote the entries of the value function matrix by Cij(t;p), i,j = 1,2. 

(i) ci t i(t;p),c 2t2 {t;p) > and Ci,i(t; p) = Ci t i(t;-p), c 2t2 (t;p) = c 2a (t;~p), sgn(ci, 2 (t; p)) = sgn(p) and 
ci,2(*;p) = -ci,2{t]-p). 

16 If n > 2, the situation is more complicated: positive correlation is not transitive in general. Hence, we cannot use Definition [XZ] 
below for the characterization of well versus poorly diversified portfolios. If the correlations of the price processes satisfy transitivity, 
some of the results can be generalized. 

17 Cross price impact and correlation can have conflicting influences on the value function and the optimal strategy. We discuss 
cross asset impact at the end of Section|4.2.2| 
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(ii) Ci.i(t;p) and c 2i2 (t;p) are strictly increasing in p on [—1,0) and strictly decreasing in p on (0, 1]. Ci t2 {t;p) 
is increasing in p on [—1,1]. 

The risk mitigati on op portunity created by a strong correlation of the price increments becomes apparent 



again in Proposition 4.7 Liquidating a single asset position x — (xi,0) T results in the cost v(t,(x,0)) 
ci.i(£; p)x\. which exhibits a strict local maximum at p — and decreases as the correlation between the two 
assets becomes stronger (irrespective of the sign of the correlation). This implies in particular that it is optimal 
for the investor to trade in both assets (unless p = 0) even if the current position in one asset is zero. 

Proposition |4 . 5| suggests that an optimal liquidation strategy never changes the sign of the asset positions of 
a well diversified portfolio and always seeks to turn a poorly diversified portfolio into a well diversified portfolio. 
The following proposition confirms this conjecture. 

Proposition 4.8. Let t G [0, T) and x g R 2 be the portfolio position at time t. 

(i) If x is well diversified, then X*(s) is well diversified for all s G [t,T) with sgn(X*(s)) = sgn(a;i) for 
i = 1,2. 

(ii) If x is poorly diversified, then sgn(Jf*(s— ) — T)*(s, X*(s— ))) 7^ sgn(xi) for s G \t,T), where 

t = inf{s > t\ sgn(X*(s)) 7^ sgn(xi) or X*(s) = for some i = 1, 2} A T > t a.s. 



By Proposition 4.8 the investor trades in both assets during the entire trading horizon [0, T] if the portfolio 
is well diversified. If the portfolio is poorly diversified the execution of the dark pool order in one of the assets 
always changes the sign of the position. If dark pool orders are never executed, it can be optimal to decrease 
the risk costs by changing the sign of the position in one of the assets by only trading at the exchange (as it is 
the case in the numerical example underlying the left picture of Figure [3]); in general, this is not the casep^j 

In the following, we specify the dependence of the optimal strategy on p. By Theorem |3.9[ we have 

£.i(t,x) = x-{c\,\{t;p)x\ +ci,2(t;p)x 2 ), &(t,x) = ^-(c2,2{t;p)x2+ci !2 (t;p)x 2 ), (14) 

Ai A2 

ril(t,x) =xi + Cl > 2 t,p x 2 , ri2(t,x) = x 2 + Cl ' 2 f ,P \ xi. (15) 

Proposition 4.9. Let t G [0, T), x = (xi, x 2 ) T , X\, x 2 > and i = 1, 2. 

(i) rji{t,x) is strictly increasing in p for p G [—1,1]. 
(ii) £j(f, x) is strictly increasing in p for p G [—1,0). 
Analog results hold for x±, x 2 < and sgn(a;i) ^ sgn^)- 



For a well diversified portfolio, the profit from diversification is increasing in \p\ (cf. Proposition 4.6 1. There- 
fore, the investor decreases her trading activity both at the exchange and in the dark pool for larger p. For a 
poorly diversified portfolio, this is not necessarily the case (cf. the left picture of Figure [2] and the corresponding 
discussion); in contrast to the optimal dark pool orders, the trader might decrease her trading intensity for large 
positive p in order to save price impact costs while waiting for the execution of an order in the dark pool (which 



yields a well diversified position by Proposition 4.8 (ii) below). We illustrate the dependence of u* on p in the 
middle and the right picture of Figure [2] In the displayed case, the trading intensity of the second asset is not 
increasing in p if x is poorly diversified (i.e., p > 0). On the other hand, the optimal dark pool order is strictly 
increasing for p G [—1,1]. The symmetry of the graph in the right picture follows directly from the symmetries 
of ci_2 and 02,2- If P = 0, the optimal strategies for the two assets are independent and follow from the formulae 
of Section |4.1[ in particular, the optimal dark pool order equals x 2 = 1. 

We deduce the general structure of the optimal strategy for initial positions Xi,x 2 > from the above 



results. If x is poorly diversified {p > 0), we have x), Tji{t, x) > by Equations (14) and (JT5J) ; hence it is 
optimal to decrease the position. This is not necessarily the case for a well diversified portfolio (p < 0); in this 
case both £i(t,x) and r]i(t,x) can be negative as ci,2 (£;/>) < (cf. also the middle picture of Figure [4] below) . 

18 Consider, e.g., a portfolio where x\ = X2, Ai = A2, tri = 02 and 9i = 62- Then the optimal trading intensities for the two 
assets must be equal until a dark pool order is executed. In particular, if the orders in the dark pool are never executed, both 
position must become zero at the same time after which further trading is not optimal. 
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Figure 3: Evolution of a portfolio consisting of two highly correlated stocks over time. The left picture illustrates 
a poorly diversified portfolio, the right picture a well diversified portfolio. In both pictures, thin lines are used 
for the less liquid first stock and thick lines for the more liquid second stock. Dashed lines correspond to trading 
without the dark pool and solid lines correspond to a realization of the liquidation process using the dark pool, 
where dark pool orders for the second stock are executed at times t\ and r 2 and for the first stock at time 
T3. Dotted lines correspond to the position which the investor aims to reach by her dark pool order for the 
respective stock. x\ = 1, X2 — 1 (left picture), x-i — —1 (right picture), T = 1, 6\ = 0.5, #2 = 3, Ai = 3, 
A 2 = 0.2, a = 4, ai = <t 2 = 1 and p = 0.9. 



For the first asset, this is the case if and only if x\ < — ci 1 (t'p) 3:2 • ^ can nence be optimal to increase the 
position in order to further reduce the risk costs of the portfolio. Note that as sgn(£j(f, x)) = sga(rji(t,x)) the 
direction of trading in the dark pool and at the primary exchange is always the same. Furthermore, it can be 
optimal to neither trade at the exchange nor in the dark pool (in the first asset) at time t if X\ = ~ ^[^w^j x 2- 
As the execution of an order in the dark pool balances the trade-off of price impact costs against risk costs, 
it should intuitively never be optimal to increase a positive position or to decrease a negative position once the 
dark pool order of one of the assets has been executed. The following proposition confirms this conjecture. 

Proposition 4.10. Let t 6 [0, T), x G K 2 be the portfolio position at time t and n be the first jump time of ir. 
Then for all s > t±, i = 1,2, 

sgn(fc(«, **(«-))) =sgn(7 7i (s,X*( S -))) =sga(X*(s-)) or X*(s-)) = m {a, X*{s-)) = 0. 

We close the section by illustrating the structure of the optimal strategy and its dependence on p by a 
numerical example. To this end, we consider two strongly positively correlated assets with Ai = 3, A2 = 0.2, 
i.e., the second asset is more liquid than the first asset. We also model the dark pool in such a way that the 
execution of orders for the second asset is more probable than for the first asset: 6\ = 0.5, 6> 2 — 3[^] We consider 
a poorly diversified portfolio x — (1, 1) T and a well diversified portfolio x = (1, — 1) T of the two stocks. 

Figure [3] shows the evolution of the two portfolios if a risk averse investor applies the optimal strategy. The 
left picture corresponds to the poorly diversified portfolio, the right one to the well diversified portfolio. In both 
cases, thin lines are used for the first stock and thick lines for the second. Dashed lines correspond to trading 
without the dark pool and the solid lines correspond to a realization of the liquidation process using the dark 
pool, where the dark pool orders for the second stock are executed at times t\ and t 2 and for the first stock 
only at time T3, i.e., dark pool orders for the more liquid stock are executed twice before any execution in the 
less liquid stock takes place. Dotted lines correspond to the position which the investor aims to reach by her 



dark pool order for the respective stock (cf. Proposition 4.3) 



For the poorly diversified portfolio, the trader tries to improve her risky position by trading out of the 



second stock (cf. Proposition 4.8 (ii) and the subsequent discussion). For this stock, trading in the primary 
venue is less expensive and being executed in the dark pool is more probable. For the well diversified portfolio, 
the portfolio position is decreasing almost linearly in time in all cases. Additionally, orders in the dark pool 
are very large for the poorly diversified portfolio and comparatively small for the well diversified portfolio, in 
line with Proposition |4.8| The reason is that dark pool orders are such that either the risk costs are decreased 
significantly by an execution (in the poorly diversified case) or they are only slightly increased (in the well 
diversified case). Note that both for the poorly and for the well diversified portfolio, these effects are stronger 
for the liquid stock; for the illiquid stock, savings in price impact costs outweigh savings in risk costs. 

19 Our choice of the parameters reflects the intuition that the asset which is more liquid at the exchange (smaller A) is also more 
liquid in the dark pool (larger 9). Theoretical findings of Ye 2011 support this choice. However, we are not aware of any empirical 
evidence for this; in some cases, the opposite parameter choice can also be plausible. 



1G 



4.2.2 Dependence on price impact 



In this section, we discuss the dependence of the value function and the optimal strategy on the price impact 
and on the cross price impact. It follows directly from the definition of the cost functional J that the value 
function is increasing in the price impact matrix A. On the other hand, the increase of v is in some sense 
bounded by the increase of A. 

Proposition 4.11. Let t £ [0,T) and x e R". 

(i) v(t,x) is increasing in A. 

(ii) Let A = diag(Aj, ; j — 1, . . . ,n). For i = 1, . . . ,n, y&*&H i s decreasing in \. 

The optimal trading strategy does not need to be monotone in the price impact parameter since risk miti- 
gation and liquidation can be conflicting desires. We illustrate this situation in Figure [4] by considering three 
different well diversified portfolios. In the first case (x\ = 0.7, x-i = —0.7; left picture) further risk mitigation 
is not optimal; the trading intensity is decreasing in the price impact. In the second case {x\ = 0.7, x-i = —1.1; 
middle picture) it is profitable to increase the position for small price impact in order to reach a position with 
even less risk costs; the optimal trading intensity is increasing and further risk mitigation is only profitable for 
small enough Ai. In the third case (x\ — 0.7, x-2 = —0.81; right picture) the optimal intensity is increasing 
for very small Ai and then decreasing; the conflict apparent in the left and the middle picture destroys the 
monotonicity in this case. 



£l(o,x) 

2.94 



0.68 




0.005 



£1(0, x) 



0.32 



0.005 



0.5 



-4.10 



£1(0, x) 



0.55 



0.36 



0.005 



Figure 4: Dependence of the optimal trading intensity for the first asset on Ai G (0.005, 0.5) for a well diversified 
portfolio x with x\ — 0.7, and different positions in the second asset: X2 = —0.7 (left picture), X2 = —1 (middle 
picture) and X2 = —0.81 (right picture). T = 1, A2 = 1, 6\ = O2 = 1, ol = 4, <t\ = 02 = 1 and p = 0.9. 

We close the section by analyzing the cross asset price impact Ai t 2 in the price impact matrix 

Ai Ai,; 



A = [ A j i 2 ) for < |Ai, 2 | < 



We first consider the case sgn(Ai.2) = sgn(p)£j We obtain the following analogs of Propositions 4.5 and 4.8 (ii). 

Proposition 4.12. Let t £ [0,T) ; sgn(Ai,2) = sgn(p) and Xi,x 2 7^ 0. 

(i) v(t, (xi,X2) T ) < v(t, (xi, — X2) T ) if and only if the portfolio x is well diversified^] 
v(t, (xi, xj)) > v(t, (x±, ~x 2 ) T ) if and only if the portfolio x is poorly diversified. 

(ii) If x is poorly diversified, then sgn(AT*(s— ) — 77* (s, X* (s— ))) 7^ sgn(xi) for s £ \t,T), where 

t = inf{s > t\ sgn(AT*(s)) 7^ sgn(xi) or X*(s) = for some i = 1,2} A T > t a.s. 



(Hi) If x is well diversified, then sgn(AT*(s— ) — i]*(s, X*(s— ))) = sgn(xi) for s G [t, t), where t is as above. 

20 This case is more intuitive than the converse: if the assets are positively correlated, buying in the first asset should rather 
increase than decrease the price of the second asset. 

), we call a portfolio x = {xi,X2) (21,2:2 7^ 0) well diversified if either the signs of the 
positions are equal (sgn(xi) = sgn(:r.2)) and p, A12 < or if the signs of the positions are different and p, A12 > 0. Otherwise, the 
portfolio is poorly diversified. 



1 Similarly as before (cf. Definition 



4.4 
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For well diversified portfolios, we can recover only a part of Proposition 4.8 (i): by (iii), it is not optimal to 
change the sign of the position by placing oversized orders in the dark pool. However, it can be optimal to turn 
a well diversified portfolio into a poorly diversified portfolio by trading in the primary exchange. We illustrate 
this in the left picture of Figure [5] in the displayed scenario, the trader holds a positive position in the second 
asset only. The trading intensity of the first asset is denoted by the solid line. If Ai^ = it is optimal to sell 
stocks in the first asset which is positively correlated to the second asset (cf. the results of Section 4.2.1 in 
particular the discussion following Proposition 4.7). If Ai,2 increases, it becomes profitable to buy stocks in the 
first asset: as the cross price impact is positive, this allows the trader to increase her trading intensity in the 
second asset (dashed line) without increasing the overall impact costs too much. In parallel, the trader places 
a sell order for the first asset in the dark pool as she does not want to sell the stocks at the exchange again 
(which would result in additional price impact costs due to the sign of Ai^)- 



£*(0,x) >7i*(0,aO ei(0,x) Ci(0,x) 




Figure 5: The left picture illustrates the dependence of the optimal trading intensity on Ai t 2 G (0, 0.7) for a 
portfolio with x\ — 0, X2 — 1. The solid line refers to the first asset, the dashed line to the second asset. 
The middle picture illustrates the dependence of the optimal dark pool order on Ai ; 2 G (—0.7, 0) for a portfolio 
x\ = 1, x<i = —1. in both pictures, T = 1, Ai = A2 = a = a\ = 02 = 1, 0\ = 82 = 3 and p — 0.2; in 
particular, sgn(p) = sgn(Ai t 2) for the left picture and sgn(p) ^ sgn(Ai ! 2) for the middle picture. The right 
picture illustrates the dependence of the optimal trading intensity on 0\ £ (0, 10) for a well diversified portfolio 
x (x\ — 0.25, X2 — —1, p — 0.9; solid line) and a poorly diversified portfolio x (x\ — x\, x-i = 1; dashed line). 
A1.2 = 0; the remaining parameters are as above. 



The case sgn(Ai.2) 7^ sgn(p) (cf. Footnote 20 1 is more complicated. In this case, it can be optimal to change 
the sign of a well diversified position (in the sense of Definition 4.4 1 by placing oversized orders in the dark pool; 
for large cross price impact, the savings in impact costs resulting from the change of the position can outweigh 
the increase risk costs. We illustrate this by a numerical example in the middle picture of Figure [5] 



4.2.3 Dependence on the execution intensities 

We conclude this section by analyzing the dependence of the value function and the optimal strategy on dark 
pool liquidity (more precisely on the intensities of the Poisson process tt) . Similarly as in the single asset case 
(cf. Proposition |4. 1 1 (i) ) , the costs are decreasing in the intensities 6*;. 

Proposition 4.13. Let t € [0, T), x S R™ and i — 1, . . . , n. Then v(t, x; 6i) is decreasing in 8i. 

The same trade-offs that can cause a non-monotone dependence of the optimal strategy on the price impact 
parameters can also give rise to a non- monotone dependence on the dark pool liquidity parameters 0i. We 
illustrate the dependence of the optimal trading intensity on 9i for a poorly (dashed line) and a well diversified 
portfolio (solid line) in the right picture of Figure [5] 



5 Proofs of the main results 
5.1 Proof of Theorem Q 

For n > 2, the second summand in the matrix differential equation C T CC is in general not linear (or quadratic), 
and ([3| is not a Riccati matrix differential equation. Furthermore, a closed form solution for the corresponding 
initial value problem is not known, and the existing theory about Riccati matrix differential equations is not 
applicable directly. 
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It turns out that appropriate upper and lower bounds for the non- linear term C T CC transform to lower 
and upper bounds (P respectively Q) for the solution of the Matrix Initial Value Problem Q and yield exis- 
tence and positive definiteness of the solution on the whole interval (— oo,T] (Theorem |3.2| . To this end, we 
require a version of a well-known comparison result for matrix Riccati differential equations, which we state in 
Appendix |Aj The main step is thus to obtain adequate matrix inequalities which enable us to transfer these 
results to the Initial Value Problem (|3]). 

For C > 0, we have < CCC. The desired upper bound of CCC is a direct consequence of the matrix 
inequality stated in the following result. 



Proposition 5.1. Let C = (0^)^=1. . M „ S lR, nx ™ be a positive definite matrix and 0i > 0, i = 1, 



,71. Then 



c<#dia g (^) =ec-\ 



(16) 



where 



Proof. We prove the inequality by induction on n. It is clear for n — 1 with equality in (16 I. Let now n > 1 
and C = (ci d )i ij=lt ... in+1 € R(™+i)x(n+i) be positive definite. Define C„ € R nxn and c € R n such that 



For z — (x, y) T G R™ x R, z ^ 0, we have 

z T Cz = x T C n x + 2x T cy + c n+ i tn+ iy 2 , z T diag (y^) 2 = ^ T diag ''" ' " 2 



'n+l 



By abuse of notation, diag ( ^f^j is used both for the diagonal nx n - matrix with . . . , in the diagonal 
and for the respective diagonal (n+l) x (n+l) - matrix. Which one we refer to is always clear from the context. 



S T (0(n + l)diag(^) -c)z 

= x T (9{n) diag {^j^j - C„)x T + 9 n+1 x T diag {^) x ~ 2x T cy + j^c n+1 ^ 



+iV ■ 



(17) 



The first summand in Equation (17 1 is non-negative by the induction hypothesis. The remainder equals 
(x - diag (—)cy) 9 n+l diag (x - — — diag (*— )cy 

0(n) T dia g( 



+ I C„+l, n +l 



7 » + l 



where the first summand is nonnegative as C n (and therefore diag ( < ^ L )) is positive definite and (9 n+1 > 0. We 
have (see, e.g., the book by Horn and Johnson 1985 , Corollary 7.7. ^ ' 



< 



9(n) 



diag 



<C~ 



(18) 



by the induction hypothesis. Moreover 

'1, 



C 



and hence 
Finally, 

finishing the proof. 



nxn C n C 

1 



C I! 



c n+l,n+l 



c T C- l c 



= T C- 1 c=^>0. 
det C_ 



(19) 



Cn+i.n+i - in—: c dia g — c ^ c„ + i.„+i - c 6„ c ^ 0, 



□ 



22 For matrices A, B with < A < B, we have < B^ 1 < A' 1 . 
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Applications of |Horn and Johnson |1985| , Corollary 7.7.4 (cf. Footnote 22 again), imply the desired bound 
for CCC. Additionally, we obtain two elementary matrix inequalities. 



Corollary 5.2. Let C > and 6i > 0, i = 1, . . . ,n. Then 

CCC < ec, 
C < ndiag (c^i) and C < tr(C)J. 



(20) 



The Matrix Inequality (20) enables us to apply Theorem A.l to the Matrix Initial Value Problem ([3| such 
that we can prove existence of a solution C of |3| on the whole interval (— oo, T] and at the same time construct 
upper and lower bounds for C via the solutions of the initial value problems in Q ; these are given explicitly by 



P(l,t)=p(l,t)I, Q(l,t) = q(l,t)I, 



where 



p(l, t) := y T + ad min coth 
q(l,t) := a/ ad max coth {^J ad m&x (T - t) + k 2 (0) 
for 9 + ad m - m > respectively ad m!iX > with 



(T-t) + Kl (Z) - 



(21) 

(22) 
(23) 



Ki(l) := arcoth ( 



+ 



v le - +ad ■ 



) > 0, k 2 (1) := arcoth ( A °"° ) 



> 



and 



for 9 = adr, 



p(l,t) := 



1 



T-t 



q{l,t) := 



1 



T-t 



(24) 



respectively ad n 



0. Note also that < p(l,t),q(l,t) < oo for all t e (— oo,T] 



Proof of Theorem \3.S\ Let C{l,t) be a solution of (|3| on some interval (t\,T]; note that there exists a local 
solution by the Picard-Lindelof theorem. The symmetry of A, E and the initial value C(l,T) = II imply that 
C(l,t) is symmetric on (ti,T]. 

Let now P := AP for P as in ([2lJ. Then P(l,t) solves P' = PK^P + 9P - ad min A, P(T) = ^A on 

(— oo, T]. As P > and P commute with A, we have P(l, t) > 0. Assume that 



{t e (ti,T] | C(l, t) is not positive definite} ^ 



(25) 



A.l 



to P := -P 



and define r := sup{t € (ti,T] \ C(l,t) is not positive definite}. As C(l,T) = II > and C(l, -) is continuous , 
there exists an e > such that C(l,t) > for t € (T — e, T] and thus t < T. We apply Theorem 
and C := -C on [r, T]. We have 

P(Z,T) = — — — A > -W = C(l,T) 



and 



P' = -PA~ 1 P + 9P + ad min A, C =-CA~ x C + CCC + aS = —CA C + 6C + (aVAZJVA + CCC - 0C) . 



Let now a; G R". Applying Corollary 5.2 to — C, we obtain 



x T (aVADVX + CCC -6C- ad min A) x = m T ( VA(£> - d mia I) VA) x + x t (CCC - 9C) (t)x > 0. 
As A > 0, Theorem 



A.l 



implies C{l,t) < P(l,t) and therefore < P{l,t) < C(l,t) on (r,T]. By continuity 
of C(l, •), we have < P(1,t) < C(/,t) and thus C(l,t) > in some neighborhood of r, a contradiction to 
Assumption (25 1. Hence, C(l,t) is positive definite on the whole interval (ii,T]. Applying Theorem A.l in the 
same way as above again, yields that we may choose t\ — — oo and that < P(l,t) < C(l,t) on (— oo,! 7 ]. A 
similar argument establishes Q as an upper bound by using < CCC instead of Inequality (20 1. □ 
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5.2 Proof of Proposition |3.4 



Before we begin, we introduce the following notation. 

Notation 5.3. (i) We denote the jump times of it by (rj)j^, where tj < Tj+i (j 6 BJ almost surely (with 
the convention tq = t). 

(ii) Given the Markovian control u*(l), the Stochastic Differential Equation (|T|) possesses a unique solution. 
We denote the process controlled by u*(l) by 

X*(l,s) := X u '^{s). 

In order to prove admissibility of u*(l), we show that -)||2 is bounded by using Gronwall's inequality 

pathwise inductively on the time-intervals [r^ A T, Tj+i A T) and interlacing the jumps (cf. Remark 2.4 1. This 



can be achieved by applying the upper and lower bounds of C(l, s) from Theorem 3.2 



Lemma 5.4. Let I > Iq for Iq as in Equation pj, t € [0, T), x € R™ be the portfolio position at time t and 
8 = ^\ 9i as before. Then the following hold. 

(i) Forse \t,T), 

X*(l,s) T C(l,s)X*(l,s) < exp(9{s~t))x T C(l,t)x a.s. 
(ii) There exists a constant K independent of I such that for all s S [t,T), ||X*(Z,s)||2 < K a.s. 
Proof. (i) Let i e IN. On {t, l < T}, X*(l, •) solves the initial value problem 

x' = -A~ 1 c(i)x, x{ n ) = x* (i, n ) 

for s G [t,;, ri+i A T). Hence, as C(l) solves the Initial Value Problem 
d 



ds 



(X*(l,s) T C(l,s)X*(l,s)) 

X*(l, s) T C(l, s)X*(l, s) + X*(l, s) T — C{l, s)X*(l, s) + X*(l, s) T C(l, s)—X*(l, s) 



ds 



ds 



ds' 



= -X*(l, s) T C{l, s)A~ 1 C{l, s)X* (I, s) + X* (I, s) T C(l, s)A~ 1 C{l, s)X*(l, s) 
+ X*{1, s) T C(l, s)C(l, s)C(l, s)X*(l, s) - aX*(l, s) T Y;X*(l, s) 
- X*{1, s) T C{l, s)A- 1 C(l, s)X*(l, s) 

< 9X*(l,s) T C(l,s)X*(l,s) 

by Corollary 5.2 and the fact that C, A , £ > 0. By Gronwall's inequality, this implies 

X*(l,s) T C(l,s)X*(l,s) < exp(d(s-T t ))x T C(l,t)x. 



(26) 



Now, on {Ti + i < T}, there exits an (almost surely) unique j = 1, . . . ,n (cf. Assumption 2.2 (ii)) such that 

X*{1, Ti+i) = X*(l, ri+i-) - T i+1 ,X*(l, n +1 -))e r 

Let r] e R. Then 

(X*(l, r l+1 -) - Vej ) T C(l, r i+1 )(X*(l, n+1 -) - nej) 



X*(l,T i+1 -) 1 C(l 7 n +1 )X*(l,T t+1 -) +r 1 2 eJC(l : T l+1 )e J - 2 V e] C(1,t 1+1 )X* (Z,r i+1 -) 



X 



"(l,T i+1 -) T C(l,T i+1 )X* (l,T i+1 -) + (j]^Cjj(l,T i+1 ) 

^ AejC(l,T l+1 )X*(l,T l+1 -)) 2 , 



y'-j ,(!■-. i) 



eJC(l,r i+1 )X*(l,T i+1 -)y 



which attains its (unique) minimum in r\ for rj = rj*(l,Ti+i,X*(l,Ti-\-i—)). Hence, 

(X*(l,Ti+i-) - ry*(Z,r i+ i,X*(Z,Ti + i-))e i ) T C(Z,T i+ i)(X*(Z,T l+ i-) - rjj(l,T i+1 ,X*(l,Ti + i-))ej) 
= rmn (X*(l,T i+1 -) - nej) T C(l,T i+1 )(X*(l,T l+1 -) - rjej) 
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and therefore 



X*(l, T l+1 ) T C(l, T i+1 )X*(l, T i+1 ) < X*(l, T, + 1 -) T C(/, T l+1 )X* (I, T i+ l-). 



Using Inequalities (26 1 and (|27[) inductively, we obtain the assertion. 



(ii) By Theorem 3.2 we have A < l/p(l, s)C(l, s) and hence, 



(27) 



X*(l,s) T AX*(l,s) < ——X*(l,s) T C(l,s)X*(l,s) 
p{l,s) 



< 



X*(l,s) T C(l,s)X*(l,s) 



(28) 



p(h,s) 
< Kx T C(l,t)x 

for a constant K independent of s by (i) and the fact that p(lo) attains its minimum in [t, T). The assertion 
follows as q{t) := lim^oo q(l,t) < oo exists (cf. Equations (23 1 and (24), cf. also Lemma [5.6| below) and 
C(l,t) < \ max q(t) (cf. Theorem pO). 

□ 



tion 



2.3 (ii) and hence Proposition 3.4 



The bound obtained in Lemma 5.4 enables us to prove that u*(l) fulfills the moment conditions of Defini- 



Proof of Proposition \3.4\ Definition 2.3 (i) and (iii) are clearly satisfied. 

Let || • || 2,2 denote the matrix norm induced by the space (lR n , || • \<i). Note that || • || 2,2 is the spectral norm 
on R nx ™ and therefore (see, e.g., |Bernstem] |2005|, Theorem 8.4.9) ||A|| 2 ,2 < ||-B|| 2 ,2 for < A < B. Using 



Theorem 3.2 and 5.4 (ii), we deduce 



IE 



U*(l,s,X*(s))\\*ds 



< E 



\*-%2\\C(l,8)\\i 2 \\X*(l,s)\\$d8 



< E 



q(l,s) 4 K 4 ds 



< OO, 



for K as in Lemma |5.4| (ii). Similarly, 

i-T 



E 



\r,*(l,s,X*(s-))\\*ds 



< E 



< E 



\\C^\\l 2 \\C(l,s)\\l4X*(l,s-)\\lds 

<nC(i,s)- 1 by Corollary |5^2] 



n 8 X s max q{l,s) s 8 

I A minP(^ S ) 8 



< OO. 



□ 



5.3 Proof of Theorem I3TB1 

We first require the following moment estimate for the controlled process. 
Lemma 5.5. Let t e [0, T), x £ R" and u £ A(t). Then 



E 



sup ||X u (s) 

t<s<T 



< 00, in particular E 



\X u (sW 2 ds 



< 00. 



Proof. Let s £ [t,T], Then by Holder's inequality and a multi-dimensional version of Jensen's inequality (see, 
e.g., Kallenberg |2002| , Lemma 3.5) we obtain 



||* tt (*)H2< INI2+ / Z(r)dr 



rj{r)dTT{r) 



<27[\\x\\$ + {s-ty / U(rmdr+ / r,(r)dn(r) 
Jt Jt 



(29) 



By Definition 2.3 (ii), it is hence sufficient to consider the last summand of Equation (29 1. To this end we 
define the compensated Poisson processes Mi(s) :— iTi(s) — 0jS, i = 1, . . . , n. We note that Poisson distributed 
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random variables have finite moments and apply Ito's isometry (note that (Mj)(s) = 0js) and Holder's inequality 
to obtain 



E 



rji(r)diii{r) 



< E 

< E 

< K. 



l dr 



E 



■E 



\r]i(r)\ 4 dr 



for a constant Ki independent of s by Definition 2.3 (ii), which finishes the proof. 



□ 



Proof of Theorem \375\ Let / > l , (t,x) € [0,T) x R™ and u = (£,77) e A(i). We apply Ito's formula (see, e.g., 
the book by|0ksendal and Sulem| |2007|) to the function w(l, t, X u (t)) = X u (t) T C(l,t)X u (t). 



rT dw 
w(l,t,x)=w(l,T,X u (T)) + I VMl,s,X«(3))Z(s) - —(l,s,X u (s))ds 



ds 



« t 1 



i=l 



<^,T,X"(T))+ y /(X«(a)^(s))dfl 

-j 1 n 



W 



2—1 

n 

( ^ 0i («;(!, s, X U ( S -)) - w(l, s, X u (s-) - m(s)ei)y S (30) 

i=i 

(i, T, X U (T)) + [ T f(X u (s),Z(s))ds +J2 f T (w(l, s, X«(s-)) - w(l, s, X u ( S -)-r ?i (s)e i ))M i (d S ), 



where Mi is the compensated Poisson process Mj(s) := 7Ti(s) — 0jS and Inequality (30) follows from Proposi- 
tion 2.5 Furthermore, we have (pathwise) equality in (30 1 if and only if u — u* A - a.s. 
Taking expectations on both sides, we obtain 

n pT 

w(l,t,x) < J(l,t,x,u) +X] E [ / s, X u (s-)) - w(l, s, X u {s-) - ^{s^ifj Mi{ds) 

i=l 

with equality if and only if u = u* P (g) A - a.s. 



(31) 



It remains to show that the stochastic integrals in Inequality (31 ) are martingales. To this end, we compute 



E 



\w(l, s, X u (s-)) - w(l, s, X u (s-) - rn(a)ei) \ 2 ds 



= E 


J 


f \2X u (s~) T C(l, s)e m (s) - misfci, 


i(«,0| 2 d« 






< E 


[J 


' 2|2X il ( S -) T Ca, S ) ei 7 ?J ( S )| 2 ds 


+ E 


r /" T 

/ 2|?7,;( s ) 2 c l , l (Z,s)| 2 ds 




< 8E 


£ \\X u (s-) T \\l\\C(l,s)\\l 2 \ Vi (s 


)| 2 d S 


+ 2E 


u 


\Vi(s)\ 4 \d,i(h 


s)\ 2 ds 


<8A 2 n 


ix ( max o(/, S ) 2 )E[ [ T \\X u (s- 


)|||ds 




r 


1 

\\v{s)\\tds 





+ 2A 2 nax ( max q(l,s) 2 )E / ||r7( S )||frfs 



(32) 



< 00 



by Definition |2.3| (ii) and Lemma 5.5 where Inequality (32) follows from Holder's inequality. As (Mj)(s) = #,;s, 
this finishes the proof. □ 
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5.4 Proof of Theorem 13.61 



We start by computing the limits of the functions p{l) and q(l) given by Equations ( 22 1 and ( 23 1 respectively ( 24 1 



Lemma 5.6. Let t G [0,T) and p{l) and q(l) as in Equations ( 22 1 and ( 23 I respectively ( 24 1 . Then 



lim ■) = p(-), lim <?(7, ■) = <?(•), 

/— >-oo (— >oo 



where p and q are given by 



p(t) 
i{t) 



T + admin coth x + ad min (T -t)J - -, 



V ad max coth (Vad max (T - i)J 
/or + admin > respectively ad max > and 

p(t) := <?(£) := 



T-£ 



(33) 
(34) 

(35) 



for 9 = ad n 



respectively ad n 



0. The convergence is compact and strictly increasing on [0,T). 



Furthermore, p(l,T),q(l,T) /* oo as I —¥ oo. 



Proof. Point- wise convergence and the formulae for the limits are straightforward by Equations ( 22 1 and ( 23 1 



respectively by Equation (24 1. Strict monotonicity follows from the fact that the initial values are strictly 
increasing in I. Finally compact convergence follows from these observations by Dini's theorem. □ 



Proof of Theorem \3.6\ (i) Note first that Theorem 3.5 implies that for fixed t 6 (— oo, T], C(l, t) is increasing 
in I on (Zo,oo) for Iq as in Equation ^ in the sense of Notation 3.1 (i). The existence of the element-wise 
limit of (C(Z,£))z>o follows directly from this monotonicity and the boundedness by Aq(t)I for q(t) as in 
Equation (p4|) respectively (35). Compact convergence follows by Dini's theorem due to the monotonicity. 



Finally, p(l,T) /* oo implies lim/_ i . 00 Cmi n (i,T) = oo. 



(ii) The inequalities in (10 1 follow directly from Lemma |5. 6 1 and (i). Furthermore, the compact convergence of 
C{l,t) (and C(l,t)) on [0, T) and the fact that C(l, ■) solves the Matrix Differential Equation (|9| implies 

C(t) = lim C(l, t) = lim f (CA^C + CCC - aS) (I, s)ds + lim C(l, 0) 



I— too 



I— too 



(CA^C + CCC - oE) (s)ds + C(0) 



and hence that C solves the Differential Equation (|9| with boundary condition lim s _j.T_ c m in(s) = oo (as 
lim-s^T- p(s) = oo). 

□ 



5.5 Proof of Theorem PI 

We start by proving the following bounds for X*(l) and X*. 
Lemma 5.7. Let t G [0, T), x G R" be the portfolio position at time t. 

(i) Let I > Iq for Iq as in Equation ^ and s G [t,T). Then (cf. Notation 5.3 (i)), 

X*{l,s) T AX*(l,s)<x T Axexp(-2[ p(l,u)du) TT 

^ Jt 1 ^ 



< x T Axexp (Q{s - if) 



T-s 



t<Ti<S 

2A maJ 

2i+ex m 



p(i,n) 



T-t 



2X U 



n 



2i+eA max J t<n<s FV ' l > 



(36) 
(37) 



where p and q are as in Equations ( 22 I and ( 23 ) respectively as in ^24 
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(ii) Let s £ [t,T). The 



X* (s) T AX* (s) < x T Ax exp (6(s -t)) ^- % TT 4^ 



U.S., 



where p and q are as in Equations ( 33 I and ( 34 ) respectively as in ( 35 I . 

Proof. We prove (i); (ii) follows by exactly the same line of reasoning with the respective bounds. Let i 6 M. 
Observe that on {n < T} 



d_ 

ds 



(X*(l,s) T AX*(l,s)) = -2X*(l,s) T C(l,s)X*(l,s) < -2p(l, s)X*(l, s) T AX*(l, s) 



for s £ [ri,Ti + i A T) by Theorem 3.2 Gronwall's inequality implies 

X*(l,s) T AX*(l,s) < X*(l,T l ) T AX*(l,T i )exp(~2 J p(l,r)drj, 

in particular 



X*{1, ( n+1 A T)-) T AX*(l, ( Ti+ i A T)—) < X*(l, n) T AX*{l, n ) exp ( 



n+iAT 



p(l, r)dr 



(38) 



This implies (cf. Inequalities (27 1 and (28) in the proof of Lemma 5.4 1 



X*(l,T l+1 ) T AX*(l 7 T l+1 ) < — -X*(l,T l+1 -) T C(l,T l+1 )X*(l,T, l+1 -) 

p{l,n+i) 

< q{ l; Tl+1 \ x* (l, r t+1 -) T AX* (I, n+ i-) a.s. 



(39) 



Using Inequalities (38) and (39 1, we obtain Inequality (|36| inductively as before. Inequality (37 1 follows from 



T-r 



2A ma: 

2i+e\„ 



which is a direct consequence of Corollary |A.2| □ 

The main step in the proof of the admissibility of u* is to show that the liquidation constraint holds 
(cf. Definition 2.3 (iv)). This is (in particular) accomplished in the following proposition. 

Proposition 5.8. Let t G [0, T) and x £ R" be the portfolio position at time t. 
(i) X*(l,-) '-Z^ X*(-) a.s. compactly on [t,T). 

(ii) I ■ \\X*{l,T)\\l a.s. and m L 1 and X*(l,T) ^ X*(T) = lim X*(s) = a.s. 

s— >T— 

Proof. (i) The spectral norm |j • || 2.2 is equivalent to the matrix maximum norm, and therefore the element-wise 
convergence results from Theorem 3.6 (i) transfer to the corresponding results for the spectral norm. 

Let t < T' < T. On {r^ < T'}, X* and X*(l) solve the respective ordinary differential equations 

X' = -?(-,X) = -A^CX, X' = -£*(!,., X) = -A^C^X 

on the interval [r», r<+i A T'). We prove that the assertion follows from the continuous dependence of 
solutions of ordinary differential equations on the right hand side and initial values. To this end, we first 
require some preliminary observations. 

For s £ [t, T'] and x, y £ R n , we have 

\\A- 1 C(s)x-A- 1 C( S )y\\ 2 < max \\Q(s)\\ 2!2 \\x - y\\ 2 =: L\\x - y\\ 2 

sG[i,T'] 
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(cf. Theorem 3.6 (ii)), i.e., for all s G [t,T'], £*(s, •) is Lipschitz continuous on It™ with Lipschitz constant 
L = L(T') independent of s. Furthermore, there exits a constant K\ > 1 such that for s £ [ijT 1 '], 
|| C(s)C(s) || 2.2 < K\. We now show by induction on i G IN that for all e > 0, there exits an ij > Zo such 
that U > k-'i and for all I > k, s E [t,n AT'], \\X*(l,a) - X*(s)|| 2 < e. 

The assertion is clear for i = 0. Let i > and e > 0. By the induction hypothesis, there exists > Zo 
such that for I > 

-L{T'-t) 

|lX*(i,r i - 1 )-X*(r i _ 1 )|| 2 <6 — . (40) 

Note that on {Tj_i > T'} the induction step is trivial. We therefore fix some ui E {t^i < T'}. Let 
now li > li-i such that for I > ij, s < n A T' (recall the uniform convergence of (C(l,s))i on [t, T'j, 
Theorem [3^6] (i)) 

-L(T'-t) 

WA-'Cil, s) A- l C(s)\\ 2 , 2 < \\A-%, 2 \\C(l, a) - C(s)\\ 2 , 2 < e ^—^ (41) 
\\C(l,s)C(l,s)-C(s)C(s)\\ 2 , 2 <^ (42) 



and 



for K as in Lemma 5.4 (ii). By the continuous dependence of solutions of systems of ordinary differential 
equations on the right hand side and initial values, we have for s E [Tj_i,Tj A T') (by Inequalities (40 1, 
and (j4TJ), 

/ p -L(T'-t) -L(T'-t) \ 

- X* Mh < {e-^- + (T' - ^^^)e^ = ^ . 



||X*(Z, fa A T')-) - X*((n A T')-)|| 2 < — . (43) 



in particular 



We can conclude by using the Inequalities ( 42 ) and (|43| : 

\\X*(l,Ti(uj) A T» - X*(n A T»|| 2 

= \\X*(l, ( Ti (u) A T')—,lS) - C{1, Ti(w) A T')C(Z, 7i(w) A T')X*(Z, (t<(w) A T')—,oj) 
- X*(( n (u;) A !")-,«) + C(n(u,) A T')C(r 4 ( W ) A T')X* (( n (u;) A I")-, w)|| a 

< ||X*(Z, (t,(w) A T')-, w) - X*((Ti(w) A 7*)-, w)|| 2 

V V ' 

<e/(3ifi)<e/3 by Inequality (43]l 

+ ||X*(J, (Ti(u) A T , )-,u)||a||C(l ) T i (w) A T')C(l,n{u) A T') - C(n(tj) A T')C{ n {u) A T')\\ 2>2 

S v ' 

<e/3 by Inequality (|42} 

+ ||(7(t 4 H A T')C(t,H A T')|| 2 ||X*^ ( Ti (w) A T')— ,w) - ^(faM A I")-, w)|| 2 

S v ' 

<e/3 by Inequality l |43"| l 

< e 
as required. 



(ii) For fixed lj E fl, we have by Lemma 5.7 (i) that 

Z.||jr(I,T, W )||»< J^^Axap^T-t))— JJ fe). (44) 



Furthermore, 



< K 9 < oo 



t<n<T ' 

for some constant if 2 independent of I; thus for almost all lo G SI, there exists a constant _fC 3 (aj) such that 

t<n{oj)<T ^ v ' IV " 
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Therefore, Inequality (44 1 implies 



Finally, Lemma 



5.7 



E[l- \\X*(l,T,u)\\l] =^0 and I- \\X*(l,T,u)\\l ^ a.s 
(ii) implies that lim ||X*(s)|| 2 = a.s., finishing the proof. 

s-T— 



□ 



We are now able to prove that u* is indeed an admissible liquidation strategy. The main step towards this 
goal is accomplished by Proposition 5.8 (ii). It remains thus to show that u* fulfills the moment conditions in 
Definition [23] (ii). 



Proof of Theorem \3.8\ Definition 2.3 (i) and (iii) are clear and (iv) follows from Proposition 5.8 (ii). 
Furthermore, we have 



E 



U(s,X*(s))\\ids 



< E 



q(s) 4 \\X*(s)\\ 4 2 ds 



< oo 



by Lemma 



5.4 



(ii) and Proposition 5.8 (i). Finally, C(s) < nC(s) 1 by Corollary 



5.2 



and thus, as 4-r admits 



a continuous extension on [t, T] , there exists a constant K independent of s such that 



\\C{ S )C{s)\\* <V||C"il JC||to < 



" 8 ALx<?(*) 8 



< K. 



Using Lemma 5.4 (ii) and Proposition 5.8 (i) again, we can deduce from Inequality (45 1 that 



E 



\ V ( S ,X*(s-)\\ 8 2 ds 



< E 



\C(s)C(s)f 2 \\X*(s-)\\lds 



< oo. 



(45) 



□ 



5.6 Proof of Theorem 13.91 

We can directly deduce compact convergence of the optimal trading intensity in the primary venue from Theo- 



rem [3^(1) and Proposition 5.8 (i). 

Corollary 5.9. Let t E [0, T) and x 6 R™ be the portfolio position at time t. Then 

-,X*(l, ■)) — >■ £*(•, X*(-)) a.s. compactly on [t,T) 

as I — > oo. 

This enables us to finally prove the main result of the article. 
Proof of Theorem\3^ We fix t £ [0, T) and x E R™. Note first that we have 

v[t, x) > lim v(l, t, x). 

I— >oo 



(46) 



For the converse inequality, let A := {tt(T) = 7r(t)} be the set of scenarios without any dark pool execution in 
[0, T] and K : (Iq, oo) x il — > [0, oo] be the following cost function: 



K(1,uj) := / (C(l,s,X*(l,s,uj)) T AC(l,s,X*(l,s,uj)) + aX*(l,s,uj) T Y,X*(l,s,uj)) ds + l\\X*(l, T, w)\ 



Then ¥[A] > and for lj E A, K(l,A) :— K(l,uS) is independent of the specific scenario uj almost surely. By 
optimality of u*(l) (Theorem |3.5| , K(l,A) is an upper bound for K(l,-) almost surely. As lim;-^ v(l, t, x) is 
bounded and F[A] > 0, there exists a constant K such that for all / > Iq, K(l,A) < K. By the dominated 
convergence theorem this implies 



lim v(l,t,x) = lim E[K{1)) =E[lim K(l)}. 

I— >oo I— foo I— too 
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By Proposition |5.8| and Corollary |5.9| the limit in the last expression exists, so Fatou's lemma yields 
lim v(l, t, x) 



> E 



= E 



lim (C (I, s, X*(l, s)) T AC {I, s, X* (I, s)) + aX* {I, s) T SI* (Z, s)) ds + lim I ■ \\X* (I, T)\\% 



(r (s, X* (s)) T Af (s, X* (s)) + aX* (s) T EI*( S )) ds 



> v(t, x). 



(47) 



The Inequalities ( |46| and (47 1 establish that u* solves the Optimization Problem ( OPT[ ) and that the value 
function is given by v. For uniqueness, let u = (£,i]),u= (£,??) € A(t,x) and /i <E (0, 1). We define the convex 
combination u = (£, ff): 

It*) = l*{s) + (1 - M)l(s), = + (1 - n)i}(s) 

for s e [t,T). Thus, I"(s) = //A u (s) + (1 - n)X a (s) and u e A(t,x). Notice that 

P ® A [u ^ u] > implies P <E> A [f ^ >0 (48) 



as else P[lim s _j.T- X u (s) ^ lim,,-^- X u (s)] > 0, a contradiction to Definition 2.3 (iv). Hence, 

l-T 



J(t,x,u)=E / f(C(r),X a (r))dr 



< E 



rf(Z(r), X u (r)) + (1 - »)f(i(r), X*{r))dr 



(49) 



/i J(t, x, u) + (1 — /i) J(i, £, w), 



where Inequality (49 I follows from the convexity of /. We have equality in Inequality (49 I if and only if u = u 

□ 



A - a.s. by strict convexity of / in the first argument and (48 1 



6 Proofs of the results of Section H] 
6.1 Proofs of the results of Section 14.11 

We first require the following elementary result. 
Lemma 6.1. Let < a < b, x > 0. Then 



d sinh(aa;) . , , sinh(a:r) 
> — . A ;, ! > (a-b)- 



dx sinh(6.x) 



sinh(6a:) 



Proof. Note first that 



d sinh(a.T) acosh(aa;) sinh(6a;) — 6cosh(6a;) sinh(ax) 



dx sinh(&.T) 



sinh (bx) 



The result follows from elementary calculus (see Kratz 201 1| for details) 



Proof of Proposition J^.l . We let t G [0, T) and compute 

A0(T - t) 



a Aecoth(§(r-t)) 

S0 CM) = 26 



4 sinh 2 (f(T-t)) 



A 

2 



< A0(cosh (| (T - t)) sinh (|(T - t)) - f (r - t) - sinh 2 (f (T - t))) < q 

20 sinh 2 (f(T-i)) 



(50) 



□ 
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for 9 > since 



cosh (-(T - t)) sinh (-(T-t))--(T-t)- sinh 2 (~(T - <)) 

= sinh (°-(T t)) ( cosh (°-(T - t)) - sinh (|(T - i))) - ^(T - t) 

l-exp(-0(T-*))\ 0, 

1 ;M --(T-i)<0 



(note that ~ (l — exp(— 2x)) — x < for x > 0). This establishes the first and the second assertion directly; the 
third assertion follows from the first equality in Equation (12 1. For the proof of (iv), we note first that 



E[X*(t; 9)} = V[ir(t) = 0] • X(t; 9) 



sinh (§(T-f)) exp (- fi) 
sinh(fT) 



(51) 



We compute for 9 > 
d 
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H[T t] sinh ( 9 -T) cosh (|( T -t))exp(-^) 



| sinh (|t) sinh (|(T - i)) exp ( - ^i) - ^ cosh (|t) sinh (°-(T - t)) exp ( - °-t 



< 



'exp ( — |i)x 
20 sinh 2 (|T) 



((T - t) sinh (|r) cosh (|(T - i)) ~ rcosh ^nh (|(T - t))) < 



We have 



by Lemma 6.1 (cf. also Equation (|50|), finishing the proof of (iv) 

E[X*(t;9) 2 } : 



, - „ sinhff (T-i)) , 
sinh(|T) 



This term is differcntiable and strictly decreasing in 9 by Lemma 6.1 (note that 9 is strictly increasing in 9). 
Thus, by Fubini's theorem, 



E 



X*{t;9) 2 dt 



0_ 



E[X*(t;6) 2 ]dt 



d_ 

09 



E[X*(t;9) 2 }dt < 0, 



establishing (v). 

Finally, we note that 



E[C (t,X*(t; 9); 9) 2 } =¥[7r(t) = 0} 



C(t;9) 2 ~^ Q . 2 _C(t;9) 2 



A 2 



-X(t;9Y 



A 2 



E[X*(t;9) 2 



by Equation (51 1. This term is differentiable and strictly decreasing in 9 as both terms are positive and strictly 
increasing in 9. Similarly as before, we deduce (vi). □ 

Proof of Proposition^-^ (i) follows directly from the cost functional J. For (ii), we compute 



< 



9 rr* ^- ld§ f cosh (°^ T - *)) sinh ~ *)) ~ s( T - 
5 1 ' J ~ 2 Ml 



<9A 



sinh 2 (f(T-t)) 



as It < 0. Monotonicity of A" follows as in the proof of Proposition 4.1 (iii). A similar calculation yields 



d 

0(aE) 



C(t;aT,) > 



as > o, finishing the proof. 



□ 
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6.2 Proofs of the results of Section 14.21 



Proof of Proposition ^ -S\ We have v(t,x — Tje,-) = (x — rjei) T C(t)(x — r/ei) which can easily be seen to be 
minimized by rj\ (t, x). □ 



Proof of Proposition ^.^ We prove the case x\,X2 > 0, p < 0, i.e., x is well diversified. Let u* be the optimal 
strategy for the initial portfolio position (xi, — x 2 ) T . For {x\, x 2 ) T , we define the strategy u £ A(t, (si, £ 2 ) T ) in 
such a way that for i = 1, 2, s > t, Xf(s) > and |X"(s)| = \Xf (s)|; this is achieved by changing the signs of 
the trading intensities and by adjusting the dark pool orders appropriately if necessary. In particular, we have 
|£t(s, X u {s))\ = \£*(s,X u (s))| and both strategies yield the same impact costs. On the other hand, the risk 
costs of u are strictly smaller as p < 0. Hence, 

v(t, (xi, -x 2 ) T ) = J(t, (xi 1 -x 2 ) T ,u*) > J(t, (xi 1 x 2 ) T 1 u) > v(t, {xi,x 2 ) T ) 

as desired. The remaining cases follow accordingly. □ 



Proof of Proposition \4-6\ We prove the case Xx,X2 > 0, p < p < and proceed similarly as in the proof of 
Proposition 4.5 Let u be the optimal strategy for p. For p, we define the strategy u £ A(t, x) in such a way 
that for i = 1,2, s > t, Xf(s) > and \X?(s)\ = \Xf(s)\. As p < p < 0, this yields 

v(t, x; p) < J(t, x, u; p) < J(t, x, u; p) = v{t, x; p). 

The remaining cases follow in the same way. □ 

For the proof of Proposition |4.7| we first require the following symmetry results of the value function. 

Lemma 6.2. Let t £ [0, T), x £ R 2 and p £ [—1,1]. Then, v(t,x;p) = v(t, — x, ; p) and v(t, (x%, X2) T ; p) = 
v(t, (xi, ~x 2 ) T ; -p). 

Proof. We have J(t,x,u) = J(t,—x,—u) and hence the first assertion follows. The second assertion follows 
from J(t, (x 1 ,x 2 ) T , (wi,w 2 ) T ;p) = J(t, (x x , -x 2 ) T , (u 1 ,-u 2 ) T ; -p). □ 

Proof of Propositio n\4-7\ (i) The first assertion follows directly from C(t) > 0. The second assertion follows 
from Lemma 16.21 as for i = 1,2, Ci^(t,p) = v(t, e^p) = v{t,e^—p) = c^^t, — p). We directly deduce the 
third assertion as 

ci,i(*; p) + ca,a(*; P) + 2cx, 2 (t, p) = v(t, (1, 1) T ; p) = v{t, (1, -1) T ; -p) = c 1A (t; p) + c 2a {t; p) - 2c xa {t; -p). 
Finally, it follows for p < (p > 0) by Proposition |4.5| that 

ci,i(t; p) + ca,a(*; p) + 2c 1>2 (t; p) = »(t, (1, 1) T ; p) < (>)v(t, (1, 1) T ; -p) = c hl (t; p) + c 2 , 2 (t; p) - 2ci, 2 (t; p) 
and therefore Ci i2 (t;p) < (ci i2 (t;p) > 0). 



(ii) For the monotonicity of Ci t i(f; •), we let p < p < and proceed similarly as in the proofs of Propositions [475] 
and 4.6 Let u be the optimal strategy for p. For p, we define the strategy u £ A(i, ei) in such a way that 
for i = 1,2, s > t, Xf{s) > and \Xf(s)\ = \Xf(s)\. As p < p < 0, this yields J(t,ei,u;p) < J(t,ei,u;p) 
with equality if and only if A 2 (s) = a.s. However, we have £ 2 (i, ei) = l/A 2 Ci j2 (i; p) < and hence 
X 2 > in some neighborhood of t with positive probability. Thus, 

ci,i(t;p) = v(t,ev,p) < J(t,et,u;p) < J(t,ex,u;p) = v(t,e\,u;p) = ci,i(t;p). 

The monotonicity of c 2j2 and the case p > follow by the same line of reasoning. 

Before we proceed, we remark that all symmetry properties and the monotonicity of ci i and c 2 . 2 in p also 
hold for C(l) (I > Iq) with exactly the same proofs. We now prove monotonicity of ci j2 (/,t;p) (I > Iq) 
first; monotonicity of Ci j2 (t;p) then follows directly from the fact that lim^oo Ci ;2 (Z, t; p) = Ci j2 (t;p). A 
straightforward computation confirms that ci, 2 (7) fulfills the following scalar initial value problem 

^, a (M;p)=^(M;p)( ^4^ + gj T^ + ^ + g 2 ) - a^p, Cl , 2 (Z,T;p) = 0. (52) 
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By the continuous differentiable dependence of C\ >2 on the parameter p, we can exchange differentiation 
with respect to t and p and obtain the following initial value problem for dC Qp '■ 

dci.2 



l ' 2 -(l,t;p) = ^^-(l,t;p)f(l,t;p) - g(l,t;p), 



dt dp 



dp 



dp 



(l,T;p)=0 



for 



, \ c hl (l,t;p) c 2 ,2(l,t;p) 
f{l,t;p) ■= . + . 



g(l,t;p) := acri<T 2 - ci,2 



l -r i7 2 , 
1 dc 1A (l,t;p) 



1 dc 2 . 2 (l,t;p) 



Op 



A? 



dp 



this implies 



dc 1}2 . 
dp 



(l,t;p) =exp(/ f(l,s;p)dsj g(l,r; p) exp f(l,s;p)dsjdr>0 



as g(l,t;p) > by (i) and the monotonicity of Cii(Z, t; •) and c 2<2 (l, t; •). 



□ 



Proof of Proposition ^ -£\ (i) We prove the assertion for 2:1,2:2 > and p < 0. Let j S IN. We assume that 
X*(jj) > (i = 1,2) on {-n, < T} (cf. Notation 5.3 1 . We compare the initial value problem for the 
controlled process on [t^Tj+i A T) (cf. Equations ( |14| )) with the case p = 0; as c\^ 2 (s; p) < 01,2(5; 0) = 
and < Ci ti (s;p) < c i:l (s;0) by Proposition Bj) X*(s;p) > X*(s;0) for s <E [tj-,Tj + i A T). For p = 
the two components of X* evolve independently according to the results of Section |4~T| In particular, the 
optimal asset position remains positive in [tj,Tj +1 AT). It follows that X*(s; p) > for all s < Tj +1 A T 
and that X*(Tj + i~) > on {tj +1 < T}. Applying Equations ( fl5| and Proposition 

{- -i • ' 

-m(T j+ uX*(r j+1 -)) = - Cl ' 2 ^ +i;p) x*( r j+1 -) > 

ci.i(T j+1 ; p) 

and the respective result for X|. The assertion now follows by induction on j. The remaining cases follow 
accordingly. 



4.7 



we obtain on 



(ii) We assume xi,x 2l p > 0. It is clear that r > t a.s. The result follows from Equations (15 1 and Proposi- 
tion [477] as for t < s < t, 

7?i( S ,X*( S )) = X*(s) + ^^4x*( S ) > X*(s). 



ci,i(s;p)' 



□ 



Proof of Proposition ^. 9\ The assertions follow directly from Equations (14 1 and (15 1 by applying Proposi 
tionHTr] 



□ 



For the proof of Proposition 4.10 we require the following results about the monotonicity of Ca and Ci 2 in 



Lemma 6.3. For i = 1,2 and t 6 [0,T), c^i(t) is increasing in t. Ci j2 (t) is decreasing in t if p > and 
increasing in t if p < 0. 

Proof. For fixed I > Iq, we consider the Initial Value Problem its solution C(l) satis fies C (l)' > Q(l)' , where 
Q{1) solves the initial value problem Q' = l/X m i n Q 2 — aa m j n I, Q{T) = II (cf. Notation 3.1 ). As in the proof of 
Proposition |3.4| Q(l) can be computed explicitly with Q(l)'{t) > for I large enough, say I > h; in particular 
d iA {l,t) > and therefore c^(i) > 0. 

For the monotonicity of ci_2 we assume p > and consider the Initial Value Problem (52) for ci, 2 (l,t). Note 
first that for all I, c[ 2 (l,T) — —apa±a 2 < 0. Furthermore, 



" (1 t\ ' n ^f ci,i(ht) , c 2i2 (l,t) \ rt / t MMi , c 2 , 2 (i,i) ^ 
c 12 (l,t) = c 12 {l,t)\^ — — 1 — \-0 1 +0 2 j +c li 2(!,t)^ — — 1 — — j 



c[,i(i,t) , a 2 ^ty 

A 2 



>0 for l>h 

This implies 2 {l,t) < for I > li and hence c[ 2 (t) < 0. The case p < follows accordingly. 



□ 
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Proof of Proposition ^. 1 0} We assume p > and that the dark pool order for the first asset is executed at time 



Ti, i.e. 



Cl,2(Tl) 

ci,i(ti) 



A^(ti); we further assume X^( T i) > 0- We have £*(ri) = and Qi T i) > 0- As 



Ci,2(s)/cj,i(s) is decreasing in s by Lemma 6.3 we have Qi s ) > ano - £i( s ) < (in particular vfeis) > and 



? yi( r i) < 0) UI1 til the next jump time of 7r. The result follows inductively as the number of jumps is almost 
surely finite. The proof for the remaining cases is analog. □ 



Proof of Proposition ^ .1 l\ The first assertion follows directly from the definition of the cost functional J. For 
the second assertion, let A, < A, . Then, 



Xi 



E 



> E 



> 



v(t, x; Xi 

aT~ 



□ 



Proof of Proposition U 



\4-12\ (i) We prove the case xi,X2 > 0, p < 0, i.e., x is well diversified, and modify the 

proof of Proposition |4.5| Let u* be the optimal strategy for the initial portfolio position (xi, —X2) T . For 
(xi,X2) T , we can define a strategy u £ A(t, (x\,X2) ) with the following properties: for s < t\ (where T\ 
is the first jump time of w) 

Ms)\ = \&'\, \X?(s)\<\Xf (s)\, (53) 
sgn(ef (s)) = sgn(£f (*)) Bgn(fi(a)) = sgn(6(*)), 
sgn(Xf («)) = sgn(X 2 u >)) => sgn(Xr( S )) = sgn(A 2 %s)). 

m needs to be defined carefully by considering all possible combinations of the signs of the trading intensi- 
ties; note that it can be necessary to change the signs of the positions. We adjust the dark pool orders in 
such a way that X"(n) £ {jX" (ti), —jXf (ti)} for some 7 < 1 (cf. the inequality in ( 53 1 ) . For s > n, 
we proceed similarly as before by defining £j(s) = 7£"*(s) (or £j(s) = — 7£"*(s)) and the dark pool orders 
as before. 

This ensures that the impact costs of u are less or equal than the impact costs of u* while the risk costs 
are strictly smaller; hence the assertion follows inductively. The remaining cases follow by the same line 
of reasoning. 



ii), (iii) As in the proof of Proposition 4.7 (i), we obtain that sgn(ci i2 (i; p, A1.2)) = sgn(Ai,2) = sgn(p). The 



assertions follow from the fact that 77* is as in Equations (15 1 (also in the presence of cross price impact). 

□ 



Proof of Proposition 4J_3 Let $4 < 9i. For any matrix C > 0, we have C{9i) > C{0i) (cf. Equation (W 
Similarly as in the proof of Theorem 3.2 we obtain that the respective solutions of the Initial Value Problem 
fulfill C(l,t;0i) < C(l,t;8i) for I large enough. The assertion follows by taking the limit / — ► 00. 



3J 
□ 



Appendix 

A Riccati matrix differential equations 

In this section, we state a well-known comparison result about matrix Riccati Equations in the form in which 
we apply it in the proof of Theorem 3.2 A standard textbook is the one by Reid 1972| . A proof for the specific 
form of the theorem can, e.g., be found in Kratz |2011| . 

Theorem A.l. Let A(t), Bp(t),Cp(t), Bq{t),CQ{i) £ H nxn be piecewise continuous on R. Furthermore, let 
B P (t), C P (t), B Q (t), C Q {t) (t £ R) and S P , S Q £ R nxn be symmetric. Let t > t x > -00 and 

Sq <S p , < B Q (-) < Bp(-), Cp(-) < C Q (-) 
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on (t2,to\. Assume that the initial value problem 

P' = —A T P — PA — PBpP + Cp, P(t ) = S P 
possesses a solution P on (ti,to\. Then the initial value problem 

Q' = -A T Q - QA - QBqQ + Cq, Q{t ) = S Q 

possesses a solution Q on (ii,io] an d P(t) ^ Qif) on (^l 5 ^o] ■ 

We apply the theorem to scalar Riccati equations with constant coefficients and obtain a useful lower bound 
for their solution. 

Corollary A. 2. Let y be the solution of the scalar initial value problem y' = y 2 + ay — b, y(T) = c, where 
a, b > 0, c > 0, b < c 2 + ac and d := a 2 /4 + b > 0. Then for t G (— oo, T], 

y(t) > 1 • 

" T-t+ -r 1 ^ 2 

c+a/2 

Proof. As d > 0, we have that the solution z of the initial value problem z' = z 2 — d, z(T) = c + a/2 fulfills 

z(t) > - - 1 1 - on {-oo,T] 

1 t+ c+a/2 



(cf. Theorem A.l compare z with the solution of /' = / , f(T) = c+ |). The assertion follows directly from 
the fact that y\J) = z{t) - a/2. ' □ 
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